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NOMENCLATURE 


Vector  Potential  function 
Psi  (40  component  of  A 

Psi  component  of  A  in  region  I 

Psi  component  of  A  in  region  II 

Psi  component  ">f  A  in  region  III 

Psi  component  of  A  in  region  IV 

Psi  component  of  a  in  region  V 

Magnetic  flux  density  or  magnetic  induction 
Magnetic  flux  density  in  medium  I 

Magnetic  flux  density  in  medium  II 

Normal  component  of  B  in  medium  1 

Normal  component  of  B  in  medium  2 

Radial  component  of  the  magnetic  flux  density 

Radial  component  of  B  in  region  I 

Radial  component  of  B  in  region  II 

Radial  component  of  B  in  region  III 

Radial  component  of  B  in  region  IV 

Radial  component  of  B  in  Region  V 

Eta  component  of  the  magnetic  flux  density 

Eta  component  of  B  in  region  I 


mi 


nin 


niv 


B 


nv 


B0 

dv 


el*e2'e3 


H 

H 


‘tl 


t2 


J 

J 


Eta  component  of  B  in  region  II 

Eta  component  of  B  in  region  III 

Eta  component  of  B  in  region  IV 

Eta  component  of  B  in  region  V 

Theta  component  of  the  magnetic  flux  density 

Elemental  volume 

Unit  normal  vector  in  radial  direction 
Unit  normal  vector  in  eta  direction 
Unit  normal  vector  in  theta  direction 
Unit  normal  vector  in  azimuthal  direction 

Metric  coefficients  for  a  prolate  spheroidal  coordinate  system 

Magnetic  field  intensity 
Externally  applied  uniform  field 

Tangential  component  of  H  in  medium  1 

Tangential  component  of  H  in  medium  2 

Magnetic  field  intensity  in  medium  1 

Magnetic  field  intensity  in  medium  2 

Electric  current 
Magnitude  of  Jg 

Electric  current  density 
Radial  component  of  J 


J 


Surface  current  density 


Eta  component  of  J 
Theta  component  of  J 
Psi  component  of  J 

Electric  current  density  of  internal  current  band 
Electric  current  density  of  external  current  band 
Magnitude  of 
Magnitude  of 

Unit  vector  normal  to  interface,  directed  from  medium  1  into 
medium  2 


Pm(cos  0) 

P 

P* 

P 

P 

Q^(cos  0) 
R.(i=l,2,3) 


r 

r,0,ip 


x,y,z 


Wn3  n4 


n,M 

M 

M 


Associated  Legendre  function  of  the  first  kind 

Variable  used  for  simplification 
Integer  from  one  to  infinity 

Associated  Legendre  function  of  the  second  kind 

Component  of  radius  vector  to  boundary  i  which  has  spherical 
symmetry 

Radius  of  spherical  coordinate  system 
Spherical  coordinates 

Distance  of  the  point  where  A  is  being  determined  from 
Rectangular  coordinates 
Constants  (specified  values  of  n) 

Prolate  spheroidal  coordinates 
Magnetic  permeability 
Relative  magnetic  permeability 


Permeability  of  free  space 


Permeability  of  medium  1 


Permeability  of  medium  2 


Variable  equal  to  cosh  n 
Variable  equal  to  cos  0 
Magnetic  scalar  potential 


Magnetic  susceptibility 


Vector  Laplacian  operator 


Psi  vector  component  of  the  vector  Laplacian  of  A  in  spherical 
or  prolate  spheroidal  coordinates 

Vector  Laplacian  of  A  in  prolate  spheroidal  coordinates 


(sinh2  n  +  sin2  0) 
2 

a  (sinh  ri  sin  0) 


-1/2 


ifi  \9ri  sinhq 


377  (sinh  n 


A,  ) 

r 


+  30  sinQ  00  (sin  0  Vj 


Vector  Laplacian  of  A  in  spherical  coordinates 


^  \  2  1  d  \  cot  0  %  \ 

3r2  r  3r  r2  302  r2  39  "  r2  sin2  0 


Scalar  Laplacian  operator 
Divergence  operator 
Curl  operator 

Radial  component  of  the  curl  of  A 


Eta  component  of  the  curl  A 


(VxA)n 

(VxA)g  Theta  component  of  the  curl  A 


EXECUTIVE  SUMMARY 


OBJECTIVE 

The  objective  of  this  theoretical  work  was  to  derive  solutions  to  static 
ferromagnetic  problems  that  include  current-carrying  coils,  uniform  inducing  fields, 
and  linear  and  homogeneous  ferromagnetic  bodies.  The  solutions  are  intended  to  be 
used  as  classical  benchmark  validation  problems  for  comparison  with  solutions  to 
ferromagnetic  problems  obtained  by  various  numerical  techniques  such  as  the  finite 
difference  method,  the  finite  element  method,  and  the  integral  equation  iterative 
solution  method. 

APPROACH 

After  deriving  the  governing  differential  equations  from  Maxwell's  equations 
for  classical  magnetostatic  field  theory,  the  method  of  separation  of  variables 
was  employed  to  obtain  the  problem  solution. 

RESULTS 

The  magnetic  induction  was  derived  for  several  configurations  of  ferromagnetic 
spherical  and  prolate  spheroidal  bodies  (hollow  and  solid)  with  internal  and/or 
external  infinitesimally  thin  spherical  and  spheroidal  current  bands,  respectively. 
The  magnetic  induction  is  presented  for  ferromagnetic  spherical  and  spheroidal 
bodies  in  a  constant  inducing  field  of  arbitrary  orientation.  The  ferromagnetic 
bodies  were  assumed  to  be  linear  and  homogeneous.  The  reduction  of  the  current 
band  problem  solutions  to  that  of  a  current  band  in  a  vacuum  is  shown  when  the 
permeability  of  the  ferromagnetic  body  is  allowed  to  approach  that  of  a  vacuum. 

The  application  of  the  superposition  principle,  to  obtain  a  total  magnetic  field 
solution  for  the  case  of  a  ferromagnetic  body  (hollow  or  solid)  surrounding  and/or 
surrounded  by  a  current  band  and  immersed  in  a  uniform  inducing  field  of  arbitrary 
direction,  is  discussed. 
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ABSTRACT 


Magnetic  induction  is  calculated  for  several  configura¬ 
tions  of  ferromagnetic  spherical  and  prolate  spheroidal 
bodies  (hollow  and  solid)  with  internal  and/or  external 
infinitesimally  thin  spherical  and  spheroidal  current  bands, 
respectively.  Magnetic  induction  is  presented  for  ferro¬ 
magnetic  spherical  and  spheroidal  bodies  in  a  constant 
inducing  field  of  arbitrary  orientation.  The  ferromagnetic 
bodies  are  assumed  to  be  linear  and  homogeneous.  The  reduc¬ 
tion  of  the  current  band  problem  solutions  to  that  of  a 
current  band  in  vacuum  is  shown  when  the  permeability  of  the 
ferromagnetic  body  is  allowed  to  approach  that  of  vacuum. 

The  application  of  the  superposition  principle  to  obtain  a 
total  magnetic  field  solution  for  the  case  of  a  ferromag¬ 
netic  body  (hollow  or  solid)  surrounding  and/or  surrounded 
by  a  current  band  and  immersed  in  a  uniform  inducing  field 
of  arbitrary  direction,  is  discussed. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  Program  Element  11221N,  Project  B0005,  Task  Area 
B0005-SL-001,  and  Work  Unit  2704-120. 

INTRODUCTION 

It  is  well  known  that  exact  analytical  solutions  of  Maxwell's  equations  using 
classical  formulation  have  been  limited  to  body  shapes  and  inhomogeneities  that 
conform  to  a  few  separable  coordinate  systems.  With  the  application  of  modern 
digital  computers  and  numerical  methods  to  obtain  solutions  of  many  magnetostatic 
field  problems  for  practical  applications,  the  need  for  classical  benchmark  valida¬ 
tion  problems  arose.  This  theoretical  report  presents  solutions  of  Maxwell's 
equations  for  magnetostatic  problems.  It  summarizes  twelve  different  problem  solu¬ 
tions  and  discusses  how  to  obtain  the  total  field  solution  to  many  others  through 
the  application  of  the  superposition  principle.  Many  of  these  problem  solutions 
may  be  used  as  benchmark  type  classical  solutions  and  for  research  in  studying 
magnetostatic  effects.  In  addition,  the  solution  techniques  and  verification 
methods  presented  in  this  report  show  the  fundamental  techniques  of  solving 
magnetostatic  boundary  value  problem  solutions  of  Laplace's  and  Poisson's  equations 
for  spherical  and  prolate  spheroidal  coordinate  systems. 
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COORDINATE  SYSTEMS 


SPHERICAL  COORDINATE  SYSTEM 

The  spherical  coordinate  system  is  formed  by  the  intersection  of  coordinate 
surfaces  of  concentric  spheres,  cones  with  apexes  at  the  center  of  the  spheres, 
and  half  planes  emerging  from  the  axis  of  the  cone.  The  three  coordinates  of  a 


point  are  the  radius  r  of  a  sphere,  the  half-angle  0  of  the  cone,  and  the  angle 
between  a  half-plane  and  the  x  axis.  Figure  1  depicts  the  spherical  coordinate  sy 
tem.  With  each  point  in  the  spherical  coordinate  system,  there  are  associated 


three  mutually  perpendicular  unit  vectors  e  ,  eft,  and  e,. 


x  =  r  sin  9  cos 
y  =  r  sin  6  sin  0 
z  =  r  cos  6 


Figure  1  -  Spherical  Coordinate  System 
and  the  Corresponding  Unit  Vectors 


■^JP 


PROLATE  SPHEROIDAL  COORDINATE  SYSTEM 

The  prolate  spheroidal  coordinate  system  can  be  formed  by  rotating  the  two- 

dimensional  elliptic  coordinate  system,  whose  traces  in  a  plane  are  confocal 

ellipses  and  hyperbolas,  about  the  major  axis  of  the  ellipse.  * 

2 

Flammer  notes  that  it  is  customary  to  make  the  z-axis  the  axis  of  revolution. 
Figure  2  depicts  the  three-dimensional  prolate  spheroidal  coordinate  system.  In 
this  case,  the  coordinate  surfaces  are:  prolate  spheroids  for  n  =  constant;  hyper¬ 
boloids  of  two  sheets  for  0  =  constant;  meridian  planes  for  tp  =  constant.  The 
prolate  spheroidal  coordinates  shown  in  Figure  2  are  related  to  rectangular 
coordinates  by  the  following  transformation  equations: 

x  =  a  sinh  n  sin  0  cos  ip  (la) 

y  =  a  sinh  n  sin  0  sin  ip  (lb) 

z  -  a  cosh  n  cos  0  (1c) 


where  0  <  n  <  °° 

0  <_  0  _<  TT 
0  _<  ip  <  2tt 

We  have  denoted  the  interfocal  distance  by  2a  and  the  prolate  spheroidal  coordinates 

by  (n,0,ip) . 
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BASIC  EQUATIONS 


FIELD  EQUATIONS 

The  formulation  of  the  present  boundary  value  problems  implies  the  solution 
of  Maxwell's  equations  for  each  medium  subject  to  the  classical  boundary  conditions. 
Starting  with  the  general  form  of  Maxwell's  equations  and  the  constitutive  relations 
between  E  and  D  and  between  B  and  H  as  given  below,  a  general  solution  may  be 
derived . 


V  x  H  = 


V  x  E  = 


9B 

at 


V  •  D  =  p 


V  •  B  =  0 


(2a)* 


D  =  e  E  +  P 
o 


B  =  p  (H  +  M) 
o 


where  H  =  magnetic  field  intensity  (A/m)** 


_  2 
J  =  electric  current  density  (A/m  ) 

—  2 
D  =  electric  flux  density  (C/m  ) 

E  =  electric  field  intensity  (V/m) 

—  2 
B  =  magnetic  flux  density  (T  or  Wb/m  ) 

3 

p  =  free  charge  density  (C/m  ) 

_  2 
P  =  polarization  (C/m  ) 

M  =  magnetization  (A/m) 

=  permittivity  of  vacuum  =  8.85  pF/m 


ViQ  =  permeability  of  vacuum  =  400tt  nH/m 


(2b) 


*The  del  operator  V  is  defined  with  respect  to  the  rectangular  coordinate  sys¬ 
tem  and  is  strictly  valid  in  a  rectangular  coordinate  system  only.  Very  often  Vx 
and  V  •  are  used  generally  as  equivalent  symbols  for  curl  and  divergence.  This  use 
is  followed  in  this  report . 

**Def init ions  of  symbols  are  given  on  page  vii. 
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For  the  magnetostatic  case,  the  applicable  Maxwell's  Equations  (2a)  reduce  to 


V  x  H  =  J  V  •  B  =  0 

and  the  constitutive  relation  from  Equations  (2b)  is 

B  =  p  (H+M) 
o 


(3a) 


(3b) 


In  general,  for  ferromagnetic  materials,  B  is  a  nonlinear  function  of  H 

B  =  f(H)  (A) 

where,  as  shown  in  Figure  3a,  B  is  not  a  single  valued  function  of  H.  The  function 
f(H)  depends  upon  the  magnetic  history  of  the  material,  that  is,  how  the  metal 
became  magnetized.  This  is  referred  to  as  hysteresis.  It  is  also  noted  that  any 
magnetic  property  of  a  ferromagnetic  material  has  meaning  only  if  it  is  considered 
together  with  its  complete  magnetic  history. 

In  certain  practical  engineering  problems,  the  variation  in  the  magnetic  field 
intensity  is  small,  and  the  functional  relationship  between  B  and  H  is  approximately 
linear  (see  Figure  3b).  For  the  linear  case  where  the  material  is  isotropic,  the 
magnetic  induction  B  is  related  to  the  field  intensity  H  by  the  relationship 


B  =  U0  (Xm+1)  H  “  Poyr  H  =  yH  (5) 

=>  magnetic  susceptibility  (dimensionless) 

=  magnetic  permeability  of  media  (H/m) 

=  =  relative  permeability  (dimensionless) 

=  permeability  of  vacuum  =  AOOtt  nH/m 

This  report  assumes  that  the  ferromagnetic  bodies  have  isotropic  and  linear  material 
properties . 


where 


*m 


(Xm+1> 
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Figure  3b  -  Curve  for  a  Ferromagnetic  Material 
at  Low  Inducing  Fields 

Figure  3  -  Typical  Magnetization  Curve 
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SUPERPOSITION  PRINCIPLE 

Maxwell's  Equations  (2a)  are  linear  partial  differential  equations.  As  a  con¬ 
sequence  of  this  linearity,  the  superposition  principle  states  that,  generally,  any 
sum  of  the  solutions  of  Maxwell's  equations  is  again  their  solution.  Combined  with 
the  uniqueness  theorem,  which  states  that  only  one  solution  of  Maxwell's  equation 
satisfies  any  set  of  prescribed  boundary  conditions,  the  superposition  principle 
justifies  any  series  or  sum  solution  of  Maxwell's  equations. 

Thus,  if  one  desires  to  find  the  magnetic  field  solution  to  a  system  consisting 
of  a  ferromagnetic  body  in  a  uniform  field  and  in  the  presence  of  current  carrying 
conductors,  the  superposit ion  principle  may  be  applied.  The  magnetic  field  solution 
for  a  ferromagnetic  body  in  a  uniform  field  only  is  obtained  first,  then  the  magne¬ 
tic  field  solution  for  the  same  ferromagnetic  body  in  the  presence  of  the  current 
carrying  conductors  only  is  determined.  The  total  magnetic  field  solution  is  then 
the  sum  of  the  two  independent  solutions.  This  technique  allows,  for  example,  one 
to  find  the  total  field  solution  for  a  hollow  prolate  spheroid  immersed  in  a  uniform 
f ield  and  surrounded  by  a  current  band . 

MAGNETIC  INDUCTION  OF  BODIES  IN  UNIFORM  FIELDS 

For  the  case  of  a  ferromagnetic  body  of  permeability  jj  in  a  uniform  field  in 
the  absence  of  current  carrying  conductors.  Maxwell's  Equations  (2a)  reduce  to 


V  x  H  =  0 


(6a) 


V  •  B  =  0 


(6b) 


Because  the  curl  of  the  gradient  of  any  scalar  function  f  is  found  to  be  identically 

zero  [7xVf=ol  the  magnetic  field  intensity  H  is  derivable  as  the  gradient  of  a 

scalar  potential  .  That  is 
m 


H  =  -V<b 

m 


(7) 
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where  <j>m  is  the  magnetic  scalar  potential  (in  amperes)  .  Using  Equations  (5)  and 
(7)  to  find  the  magnetic  flux  density  and  substituting  the  result  into  Equation  (6b) 
reduces  to 

V2cf>  =  0  (8) 

m 


which  is  known  as  Laplace's  equation.  This  is  the  governing  differential  equation 
for  the  problem  of  a  body  immersed  in  a  uniform  field. 

The  general  boundary  conditions  to  be  satisfied  at  the  interface  of  dissimilar 
materials  may  be  derived  from  the  limiting  integral  form  of  Maxwell's  equations  and 
are  given  by 


n  •  fB  -B  ) 
12  ^2  l’ 


or  B 


nl 


9<t> 


Bn2  °r 


ml 


9<J> 


m2 


1  3n 


2  9n 


(9a) 


n12  x  (H2~H^) 


=  0 


or  Ht2  =  Htl  or  4>ml  =  ^ 


(9b) 


where  the  subscripts  1  and  2  indicate  the  media  under  consideration,  and  n^2  denotes 
the  unit  vector  normal  to  the  interface  and  directed  from  medium  1  into  medium  2. 

MAGNETIC  INDUCTION  OF  BODIES  DUE  TO  CURRENT  CARRYING  CONDUCTORS 

The  divergenceless  (V.B=0)  nature  of  the  magnetic  flux  density  in  conjunction 
with  the  fact  that  the  divergence  of  the  curl  of  any  vector  function  is  zero 
[V’(VxF)=o]  allows  the  introduction  of  the  magnetic  vector  potential  field  (A) 


B  =  VxA 


(10) 
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■pp 


where  A  is  the  magnetostatic  vector  potential  function  in  webers  per  meter.  The 
substitution  of  Equation  (10)  into  Equation  (3a)  gives  the  fundamental  equation  of 
the  vector  potential  of  the  magnetostatic  field. 


1 

y 


v  x 


(VxA)  -  (VxA)  xV 


_1 

y 


j 


(ID 


For  homogeneous  materials,  as  assumed  in  this  report,  the  magnetic  permeability  is 
spatially  invariant.  Hence 


(12) 


and  Equation  (11)  reduces  to 


V  x  V  x  A  =  MJ 


(13) 


Using  the  vector  identity 


V  x  V  x  A  =  V  (V-A)  -  $  A 


(14) 


Equation  (13)  becomes 


V(V-A)  -  4*  A  =  MJ  (15) 

The  magnetostatic  vector  potential  is  characterized  by  the  important  property  that 
its  divergence  can  be  conveniently  chosen  to  be  zero. 

V-A  =  0  (16) 
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Equation  (15)  reduces  to  the  vector  Poisson  differential  equation. 


$  A  =  -  yj  (17) 

This  is  the  governing  equation  for  our  calculations. 

The  general  boundary  conditions  to  be  satisfied  at  the  interfaces  of  stationary 
dissimilar  media  may  be  derived  from  the  limiting  integral  forms  of  Maxwell's  equa¬ 
tions  and  are  given  by 

n12  *  (B2-Bl)  =  0  or  Bnl  =  Bn2  (18a) 

“12  *  (VV  "  '  Js  (18b> 

where  the  subscripts  1  and  2  indicate  the  media  under  consideration,  and  n^  denotes 
the  unit  vector  normal  to  the  interface  and  is  directed  from  medium  1  into  medium  2. 
In  the  case  where  the  materials  are  linear  and  isotropic,  Equations  (18a)  and  (18b) 
become 


n12  *  (ll2H2_lJlHl)  0 


B2  B1 

ni2  x  l  v2  yx J  Js 


(18c) 


(18d) 


where  Jg  is  a  true  surface  current  density  that  may  exist  at  the  interface.  At  an 
interface  where  Jg  =  0,  Equations  (18b)  and  (18d)  need  to  be  modified  accordingly. 
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SOLUTIONS  FOR  SPHERICAL  BODIES 


SOLID  SPHERE  OR  SPHERICAL  SHELL  IN  A  UNIFORM  FIELD  OF  ARBITRARY  DIRECTION 

Several  important  types  of  problems  relating  to  magnetized  bodies  in  an  exter- 

3  4  5  i 

nal  magnetic  field  have  been  solved  ’  ’  by  determining  the  solution  to  Laplace  s 

equation  for  the  magnetic  scalar  potential.  Generally,  these  solutions  have  been 

derived  for  the  case  of  the  uniform  external  magnetic  field  in  the  direction  of  the 

z  axis  of  a  spherical  coordinate  system.  Both  constant  external  field  problems, 

solid  and  shell,  were  solved  and  programmed  on  the  digital  computer  by  D.A.  Nixon  of 

the  Center,  for  the  case  of  an  arbitrarily  orientated  external  magnetic  field. ^ 

The  solutions  found  in  Reference  6  were  presented  in  Cartesian  coordinates.  The 

problem  of  finding  the  magnetic  induction  for  an  infinitesimally  thin  current  band 

surrounding  a  spherical  shell^  can  be  generalized  to  include  an  external  magnetic 

field.  Linear  superposition  may  be  applied  to  find  the  solution  in  this  case. 

Therefore,  in  Appendix  A,  the  Cartesian  expressions  were  converted  to  spherical 

coordinates  to  be  compatible  with  other  problem  solutions  in  this  section  of  the 

report . 

SOLID  SPHERE  SURROUNDED  BY  AN  INFINITESIMALLY  THIN  SPHERICAL  CURRENT  BAND 

We  now  solve  the  boundary  problem  of  a  ferromagnetic  sphere  of  radius  R^  and 
homogeneous  permeability  surrounded  by  an  infinitesimally  thin  current  band  of 
radius  R£  having  a  constant  current  density  J.  Figure  4  identifies  the  three 
regions  of  interest.  Regions  II  and  III  have  a  permeability  equal  to  the 
permeability  of  vacuum  y  which  for  convenience  will  be  labeled  y^.  The  problem's 
spherical  symmetry  suggests  that  a  spherical  coordinate  system  such  as  that  shown 
in  Figure  1  be  used  in  the  problem  solution. 
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Ampere's  law  states 


V  x  H  =  J  (19) 

and,  because  V  •  B  =  0,  the  induction  B  must  be  the  curl  of  some  vector  field  A. 

The  governing  differential  equation  for  A,  when  homogeneous  and  linear  materials  are 
considered,  is  from  Equation  (11). 


* 

$  A  =  -pj 


(20) 


-  ■  —  —  —  n 

*We  note  that  a  distinction  is  drawn  between  the  operator  called  the  scalar 
Laplacian  operator  and  the  vector  Laplacian  operator  designated  by  0  .  The  vector 

Poisson's  equation  in  rectangular  coordinates  can  be  treated  as  three  uncoupled 
scalar  equations  as  shown  below. 


A 

6 

32A  32A  32A 

X  X  x 

A 

+  e 

32Ay  32Ay  82Ay 

X 

_3x2  3y2  3z2  _ 

y 

_3x2  3y2  9z2  _ 

4-  P 


32A 


3x 


32A 


3y 


32A 


3z 


=  e  J  + 

X  X 


<*  J  +  e  J 
y  y  z  z 


where  $A  =  V^A^  =  for  i  =  x,y,z.  However,  if  the  vector  Poisson's  equation 
is  resolved  into  orthogonal  components  in  other  coordinate  systems,  the  differential 
operation  mixes  the  components  together  giving  coupled  equations  as  shown  below  for 
spherical  coordinates. 


1A 


32A  ,  3A  ,  92A  q  9A 

r  2  r  21  r,  cot  0  r 

A  '  r  .  2  r  3r  ~  2  Ar  2  .  2  2  30 

9r  r  r  30  r 


,  3  A  0  3A0  .  Q  0  3A. 

1  r  2  0  2  cot  6  2  ip 

2  ,  2  Q  ,,2  ‘  2  36  '  2  0  “  2  .  Q  3^ 

r  sin  0  3v  r  r  r  sin  0 


2  3A0 

A0 

+  1  3V 

cot  0  ^A0 

i  3% 

a  2 

l_3r 

r  3r 

2 

r 

sin2 

+  2  2  + 

6  r  30^ 

2  30  2 

r  r 

sin2  0  3^ 

3A 

r 

2  cot  0 

% 

i  -s* 

32A.  „  3A, 

_ ±  ,  2  _J i 

1  » 

+  i  3\ 

39 

2  ‘  ft 

r  sin  0 

3ip 

+ 

~  2  r  3r 

L3r 

2  •  2  o  v 

r  s  m  0 

+  2  2 
r  36 

+  cot_!%+ _ 1 _  Mfe  + _ 2 _ ^  +  jLcnt  L. 

2  30  2  .  2  Q  ~  2  .  Q  dip  2  .  Q  dip 

r  r  sin  0  3^  r  sm  0  r  sin  9  r 


=  J  e  +  JQe„+  J ,  e , 

r  r  0  0  ip  ip 


The  general  expression  in  spherical  coordinates  for  a  current  density  is 


J  =  e  J  +  e.J  +  e,J 

r  r  0  0  ip  ip 


f 


where  the  e  defines  unit  orthogonal  vectors.  For  stationary  currents  in  vacuum, 
the  vector  potential  function  that  satisfies  Equation  (20)  is  given  by 


A  = 


dv 


(22) 


where  dv  =  elemental  volume  in  the  current-carrying  region 

r '  =  distance  between  the  field  point  where  A  is  being  determined  and  dv  at 
the  source  point. 

From  Equation  (22)  we  see  that  the  elemental  vector  potential  dA,  due  to  a  cur¬ 
rent  element  Jdv,  is  in  the  same  direction  as  J.  It  is  well  known  from  this  fact 
that  the  lines  of  the  magnetic  vector  potential  A  are  circles  centered  about  the 
coil  or  loop  axis.  The  magnitude  of  A  along  such  a  circle  is  constant,  which  means 
that  A  is  a  function  of  the  spherical  coordinates  r  and  0  only.  Therefore,  we  know 
in  advance  for  this  problem  that  A^  is  the  only  component  of  A  existing  at  the  field 
point.  The  infinitesimally  thin  band  of  current,  shown  in  Figure  4,  has  only  an 
azimuthal  or  component,  which  is  a  function  of  r  and  0,  and  lies  on  the  boundary 
between  regions  II  and  III  (i.e.,  r=R2)  .  For  this  current.  Equation  (21)  reduces  to 


J  = 


e, J  (0) 
■P  'P 


if  0  <0  or  6  >e2 
if  ei  1  9  1  02 


(23) 


Therefore,  Equation  (20)  has  only  an  azimuthal  component  and  can  be  expressed  as 


$A  =  4?A,(r,0)  =  0  (in  regions  I  through  III) 

r  V 


(24) 
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When  the  vector  Laplacian  XjX  is  expanded  in  spherical  coordinates.  Equation  (24) 
can  be  written  as 


d  \  2  ^  \  cot  6  ^  _  in  regions  I 

„  2  +  r  ~  2„q2+  2  36  ~  2  .  2  Q  ~  °  through  III 

3r  dr  r  39  r  r  sin  6 


(25) 


In  order  to  solve  Equation  (25)  it  is  necessary  to  obtain  the  general  solution  in 

2 

regions  I  through  III.  Thus,  by  multiplying  Equation  (25)  by  r  we  obtain 


2  2 
r  3A 


,  2r3A, 

A  + 


3r 


3r 


3\  „  sAt 

7T-  +  9  M 

ou 


•  2  A 
sin  6 


=  0 


(26) 


Applying  the  method  of  separation  of  variables,  let  us  assume  that  A^ 
expressed  as  a  product  of  two  functions 


can  be 


A^  =  R(r)0(6) 


(27) 


where  R(r)  is  a  function  of  r  only  and  0(0)  of  0  only.  Substituting  this  form  of 
the  vector  potential  A^  into  Equation  (26),  we  have,  after  separation  of  variables 


d2R(r)  2  dR(r )  _  p(p+l)R(r) 

,2  r  dr  2 

dr  r 


(28a) 
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+ c»c  e  ♦ 


d0 


de 


p(p+i)  - 


sin^  0 


0(6)  =  0 


(28b) 


where  the  separation  constant  is  p(p+l)  and  p  is  an  integer  from  one  to  infinity. 
The  differential  equation 


d2Q  .A  dQ  . 
^+cot6^  + 


p  (p+1)  - 


sin^  9 


0=0 


(29) 


has,  as  a  general  solution. 


0(6)  =  0(0) 


CLPm  (cos  0)  +  D  Qm  (cos  0) 
P  p  P> 


(30) 


2 

Comparison  of  Equations  (28b)  and  (29)  shows  that  in  Equation  (29)  m  =  1.  This 
requires  that  m  always  be  unity.  The  solutions  of  Equations  (28a)  and  (28b)  are 
then  expressed  as 


R(r)  =  R  (r)  =  A  rP 
P  P 


+  B'r"(P+]) 

P 


Ola) 


0(0)  =  0  (0)  =  C  P1 (cos  0)  +  D  Q1 (cos  0) 
P  P  P  P  P 


(31b) 


The  associated  Legendre  functions  of  the  first  and  second  kind  are  designated  as 
pm(cos  0)  and  Qm  (cos  0),  respectively.  Therefore,  the  general  solution  of 
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Equation  (25)  in  regions  I  through  III  may  be  formed  from  the  product  of  the  solu¬ 
tions  in  Equation  (31) ,  which  yields 


A  =  R(r)0(9)  =  £  R  (r)0  (Q)  (32a) 

V  P  P 

p=i 


oo 

E 

P=1 


(c  p1 (cos  e)+D  Q1 (cos  0) 
\  P  P  P  P 


(32b) 


In  the  spherical  case,  associated  Legendre  functions  of  the  second  kind  are  infinite 
at  cos  9  =  ±1,  and,  thus,  cannot  be  included  when  the  region  under  consideration 
includes  the  symmetry  axis.  Therefore,  the  constant  must  be  set  equal  to  zero, 
and  Equation  (32)  reduces  to 


P^(cos  0) 
P 


(33) 


where  A  =  A  C  and  B  =  b'  C  . 

P  P  P  P  P  P 

The  form  of  the  potential  in  each  of  the  regions  (I  through  III)  is  determined 
from  Equation  (33).  These  magnetostatic  vector  potentials  in  regions  I  through  III 
are: 


£  (V^)  PP(C0S  9) 


(34a) 
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=  Z  A  0rP+  — r^-iT  P1(cos  0) 


Att  =  A,tt  =  2_j  A  0r*+  P 

II  iJjII  p2  r(p+l)  P 


“  B  ^  , 

A  =  A  =  V*  ~  ,p  .  P  (cos  6) 
III  ij)III  ^  r(p+l)  p 


(34b) 


(34c) 


where,  for  the  A^  component,  =  0  because  at  r  =  0  the  potential  must  be  finite 

and,  for  the  A^^  component,  A^  =  0  because  as  r  approaches  infinity  the  potential 
must  remain  finite. 

At  each  interface,  the  basic  laws  of  magnetostatics  in  Equations  (3a)  reduce  to 
boundary  conditions  on  B  and  H  that  can  be  used  to  evaluate  the  four  constants  in 
Equation  (34).  From  Equations  (18a),  the  normal  component  of  B  across  each  boundary 
must  be  continuous,  i.e.,  (82-B^  ) *  n12  =  0  where  the  quantity  n^2  is  the  unit  outward 
normal  to  the  surface.  This  provides  the  following  boundary  conditions  which  must 
be  satisfied  by  the  solution  in  Equation  (34)  for  each  region. 


Brl  =  Brll  at  r  =  Rl 


(35a) 


BrII  =  BrIII  at  r  =  R2 


(35b) 


The  normal  component  of  the  magnetic  field  B^  is  expressed  in  terms  of  the 
vector  potential  as 


B  =  (VxA) 
r  r 


(36a) 
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(36b) 


B 

r 


r  sin  0 


—  (sin  0  A  ) 

96  V 


e 

r 


e0r 


e , r  sin  0 


where 


B  =  VxA  =  —  — - 

r  sin  0 


9 

9r 


9 

90 


9_ 

9\Jj 


0  0 


A.r  sin  0 


However,  because  the  vector  potentials  in  each  region  are  functions  of  (cos 
we  can  simplify  Equation  (35)  to  constraints  on  A, 


II 


at  r 


II 


AIII  at  r  =  R2 


The  second  set  of  boundary  conditions  is  obtained  from  Equation  (18b).  The  tan 
gential  component  of  H  across  each  boundary  must  satisfy  the  relationship 


n12  X  (H2-IV  ~  J s 


0) 


(37a) 


(37b) 


(38) 
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V-- 


where  .1  (which  equals  J(t3))  is  the  real  surface  current  density  in  the  limit  of 
s 

vanishing  width  between  the  two  regions.  Using  the  relationship  B  =  pH,  Equation 
(38)  may  be  expressed  as 


01 


=  -K6) 


(39) 


Referring  to  the  curl  in  Equation  (36),  we  can  write  Bft  as 


•  <?XA>0  ■  -f^K] 


(AO) 


From  Equations  (38),  (39),  ana  (AO),  the  tangential  components  in  regions  I 
through  III  must  satisfy  the  relationships 


—  —  4-  (rA  )  +  —  —  I—  (rA  )  =  0  at  r  =  R1 
r  3r  II  u2  r  3r  I  1 


(Ala) 


7  h  wm> +  77  7  I?  (rAn>  ■ J(e>  “  r  -  r2 


(Alb) 


Tlu-  general  expressions  for  the  potentials  in  each  region  (Equation  (3A))  are 
then  substituted  into  the  boundary  conditions  (Equations  (37)  and  ( A 1 ) )  and  solved 
ter  the  constants  A^  and  B  ^ .  There  are  four  algebraic  equations  with  four 
unknowns  and  the  potential  in  each  region  can  then  be  specifically  determined.  The 
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four  boundary  value  equations  that  must  be  solved  for  the  coefficients  are  given 

below  (where  the  index  p  is  odd  only  and  understood  to  take  on  values  from  1  to  ®). 

It  is  noted  that  the  current  J^(0)  must  be  expanded  into  a  set  of  associated 

Legendre  functions  in  order  to  evaluate  the  constants  A  .  and  B  . .  The  detailed 

pi  Pi 

expansion  is  in  the  section  of  Reference  7  entitled  "Expansion  of  the  Current 
(J^(0))  in  Associated  Lengendre  Polynomials." 


plftl 


VR1  +Bp2Rl 


-(p+D 


(42a) 


+BP2R2 


-(p+1) 


(42b) 


L 

Pi 


(p-1}  -  -  _PR  -(pf2) 
p2  1 


Ap2(p+1)R1'1'  *'  -  B^,pR. 


+ 


APi(p+1)R: 


(p-D 


=  o 


(42c) 


k[ 


~(p+2)i  .  i  r.  (p-D 


-(P+2)l  =  _ E 


BP3pR2  J  +  ^  [Ap2(P+1)R2VF^/  -  PBP2R2  j  „1 


J  (6) 


P  (cos  0) 
P 


(42d) 


The  solution  of  these  equations  to  obtain  the  coefficients  yields  (a  detailed 
derivation  is  given  in  Appendix  B) : 
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The  coefficients  A  .  and  B  ,  can  be  determined  from  Equations  (43),  and  Equations 
pi  pi 

(34)  can  now  be  used  to  specify  the  potentials  A^,  A^,  and  A^j^  in  regions  I  through 
III.  The  normal  (B^)  and  tangential  (Bq)  components  of  the  magnetic  induction  in 
regions  I  through  III  can  be  determined,  by  using  Equations  (36)  and  (40),  to  be: 
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in  the  outer  region. 


The  magnetic  vector  potential,  ^  in  the  inner  region  and  A^^ 
are  derived  for  an  infinitesimally  thin  current  band  in  a  homogeneous  medium  of 
permeability  u, ,  in  Reference  7.  The  coefficients  A  .  (i  =  1,2)  and  B  . 

r  1  1  ’  pi  ’  pi 

(i  =  2,3),  for  the  vector  potentials  for  the  present  ferromagnetic  sphere  problem, 
reduce  to  the  coefficients  of  the  potentials  in  the  two  regions  of  the  simple  current 
band  problem  when  the  permeability  of  the  sphere  p9  approaches  that  of  the  surround¬ 
ing  medium  p..  This  shows  that  the  solutions  of  the  above  ferromagnetic  current 
problem  have  the  correct  mathematical  form. 

We  note  that  derivation  of  the  solution  for  the  problem  of  a  ferromagnetic 
sphere  surrounded  by  a  coil  of  finite  width  is  found  in  Reference  7,  See  Figure  5 
for  the  geometry  of  the  problem.  The  magnetic  induction  for  this  case  is: 
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Figure  5  -  Ferromagnetic  Sphere  Surrounded  by  a 
Coil  of  Finite  Width  (yz  plane) 
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The  mathematical  solution  for  B^  in  terms  of  known  quantities,  obtained  in 
Appendix  E  of  Reference  7,  is  given  by 
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The  numerical  values  for  the  other  coefficients  can  be  obtained  from  the  equa¬ 


tions 


(47a) 


[Z]  +  K 


t  t 

P 


v  ■  M  V 


p2 


B  0R-, 

p2  1 


-(2p+l) 


(47b) 
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B  -  A  ,r,Cp+1>  +  b  ,  +  -!-2 - E 

p4  p3  3  p3  (P-2)  (P+3) 


(47e) 


SPHERICAL  SHELL  SURROUNDING  AN  INFINITES IMALLY  THIN  SPHERICAL  CURRENT  BAND 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic  spherical 
shell  of  outer  radius  R^,  inner  radius  R^,  of  homogeneous  permeability  surround¬ 
ing  an  infinitesimally  thin  current  band  of  radius  R^  and  having  a  current 
density  J.  A  constant  density  J  is  assumed.  Figure  6  shows  the  four  regions  of 
interest.  Regions  I,  II,  and  IV  have  a  permeability  equal  to  the  permeability  of 
vacuum  Mq  which,  for  convenience,  will  be  labeled  p^.  The  problem's  spherical 
symmetry  suggest  that  a  spherical  coordinate  system,  such  as  that  shown  in  Figure  1, 
be  used  in  the  solution. 

The  details  of  treating  problems  of  this  type  with  spherical  symmetry  are  dis¬ 
cussed  in  the  previous  section  of  this  report  and  in  Reference  7. 
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Figure  6  -  Infinitesimally  Thin  Current  Band  Surrounded 
by  a  Ferromagnetic  Spherical  Shell 
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The  partial  differential  equation  that  governs  this  problem  is  the  azimuthal 
component  of  the  vector  Laplace's  equation. 


4*  =  4*  A  (r,t3)  =  0  in  regions  I  through 


When  the  vector  Laplacian  is  expanded  in  spherical  coordinates,  Equation  (48)  can  be 
written  as 


3^A,  r,  3A.  3^A  3a,  A. 

2  (ft  1  $  cot  6  (ft  (ft 

2  r  3r  2  „02  2  30  ~  2  .  2 

9r  r  30  r  r  sm 


=  0  (in  regions  I 
through  IV) 


The  general  solution  of  the  equation  has  the  form 


A.  *  22  lA  rP  +  P  )  P  (cos  0) 

*  ,p.l\p  r (p+1)/  p 


where  A  and  B  are  constants  and  P  (cos  0)  is  the  associated  Legendre  function  of 
P  P  P 

the  first  kind.  The  magnetostatic  components  of  the  vector  potential  in  regions  I 
through  IV  are 
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v  •  £  (vpK (cos  0) 
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W  ^  (vr  +  pp  (cos  0) 


(51b) 


32 


AiJjIII 
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r<P+1)/ 


(cos  6) 


(51c) 
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OO 


(cos  9) 
P 


(51d) 


where,  for  the  A^  component,  B  ^  =  0,  because  at  r  =  0  the  potential  must  be  finite, 
and,  for  the  component,  A  ^  =  0,  because  as  r  approaches  infinity  the  potential 

must  remain  finite. 

At  each  interface,  the  basic  laws  of  magnetostatics  reduce  to  boundary  condi¬ 
tions  on  B  and  H  that  can  be  used  to  evaluate  the  six  constants  in  Equations  (51). 
The  first  boundary  condition  states  that  the  normal  component  of  B  across  each 
boundary  must  be  equal  to  (B^B^  •  n12  =  0.  The  vector  quantity  n12  is  the  unit 
outward  normal  to  the  surface  (in  the  spherical  case  ef) .  Thus,  the  following 
boundary  conditions  must  be  satisfied  by  the  four  regions  of  the  ferromagnetic 
spherical  shell  problem. 


B  T  =  B  XT  at  r  =  R, 

rl  rll  1 


(52a) 


Brll  =  \m  at  r  =  R2 


(52b) 


B  tti  =  B  tv  at  r  =  R, 

rill  rIV  i 


(52c) 


The  second  boundary  condition  states  that  the  tangential  component  of  H  across  each 
boundary  must  satisfy  the  relationship 


ni2  x  ~  Js 


(53) 
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where  Jg  (which  equals  J(0)  in  our  case)  is  the  true  surface  current  density  in  the 
limit  of  the  vanishing  width  between  the  two  regions.  Using  the  linear  relationship 
B  =  pH,  Equation  (53)  can  be  expressed  in  spherical  coordinates  as: 


J(6) 


(54) 


The  general  expressions  for  the  components  of  the  vector  potentials  in  each 
region  (Equations  (51))  are  then  substituted  into  the  boundary  conditions  (Equa¬ 
tions  (52)  and  (54))  and  solved  for  the  constants  A  ^  and  B^.  There  are  six 
algebraic  equations  with  six  unknowns,  thus  enabling  the  potential  in  each  region 
to  be  specifically  determined.  The  six  boundary  value  equations  that  must  be  solved 
for  the  coefficients  are  given  below  (where  the  index  p  is  odd  only  and  understood 
to  take  on  values  from  1  to  “). 

As  in  the  previous  section,  the  component  of  the  current  J^(0)  must  be  expanded 
into  a  set  of  associated  Legendre  functions  in  order  to  evaluate  the  constants  A^ 
and  B  (The  detailed  expansion  of  the  azimuthal  component  of  the  current  density 

is  given  in  Reference  7) . 


A 

P1 


_P 

“1 


+  B 


P2 


-(p+1)' 

1 


(55a) 


R2  +  VR2~<P+1>]  ‘  [Ap3R2  + 

[V?  +  Bp3R3<P+1)]  *  V[R3<P+1)] 


] 


(55b) 


(55c) 


34 


-k[ 


Ap2(p+1)RIP'1)  -  pBp28 


*  4tP 


v<p+i,Rr1 


Jp(6) 


p3(cos  0) 
P 


(55d) 


U' 


-  '"Sp382<P+2>]  +^[V<P+1)82P'1)  -  pBp2R2'(P+2)] 


(55e) 


—  ~pB  ,r:(p+2)  +  —  A  ,,(p+l)R, 

L  P4  3  J  P2  L  p3 


(p-1) 

3 


-  PBp383(Pf2>]  ’  0 


(55f) 


OO 

["Note:  JT ,  (0 )  =  J  V  K  P  (cos  0)  =  £  J  (0)  where  K  =  0  for  p  even."] 
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p=l 

The  coefficient  A^  in  terms  of  known  quantities  is  expressed  as: 
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(56c) 
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The  numerical  values  for  the  other  five  coefficients  can  be  obtained  from  the 
following  equations 


(57a) 
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The  coefficients  A  .  and  B  .  can  be  determined  from  Equations  (56)  and  (57).  Equa- 
p  i  pi 

tions  (51)  can  now  be  used  to  specify  the  potentials  A^,  Ajj-,  AIII’  an^  AIV 
regions  I  through  IV.  The  normal  (B^)  and  tangential  (B^)  component  of  the  magnetic 
induction  in  regions  I  through  IV  can  be  determined,  by  using  Equations  (36)  and 


(40)  as: 
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BrIV  =  S  Tin—  [Bp4r'(P+2>]  le  [Sln  9  PJ  (C°S  9)]  08h) 

P‘1 

The  coefficients  (i  =  1,2,3)  and  B  ^  (i  =  2,3,4),  for  the  vector  potential  of 
the  present  ferromagnetic  shell  problem,  reduce  to  the  coefficients  of  the  potentials 
in  the  two  regions  of  a  simple  current  band  when  the  permeability  of  the  shell 
approaches  that  of  the  surrounding  medium  p^  (see  Appendix  C) . 

Lebedev  et  al . ,  present  the  magnetic  vector  potential  due  to  a  dc  current  1 
flowing  in  a  filamentary  circular  loop  of  radius  rQ  inside  a  hollow  spherical  shell 
made  from  material  of  magnetic  permeability  p  (in  Reference  3  see  "The  Fourier 
Method,"  page  99).  The  components  of  the  magnetic  vector  potential  were  given  as: 
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where  r  >  R9 ;  see  Figure  7. 

c  =  2.998  x  10B  m/sec. 


Note:  Lebedev's  equations  are  expressed  in  Gaussian  units. 
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Figure  7  -  Ferromagnetic  Spherical  Shell  Surrounding 
a  Filimentary  Current  Loop 


Although  we  did  not  derive  this  solution,  a  solution  for  this  type  of  problem 
could  be  obtained  by  allowing  the  infinitesimally  thin  current  band  in  the  preeodin 
problem  to  degenerate  to  a  filamentary  current  loop  as  was  done  in  Appendix  B  of 
Reference  7 , 

SPHERICAL  SHELL  SURROUNDED  BY  AN  INFINITESIMALLY  THIN  SPHERICAL  CURRENT  BAND 

The  boundary  value  problem  of  a  ferromagnetic  spherical  shell  of  outer  radio.-; 
R., ,  inner  radius  R^ ,  and  a  homogeneous  permeability  p2>  surrounded  by  an  iniinitesj 
mally  thin  current  band  of  radius  R^  having  a  current  density  J  was  solved  in 
Reference  7.  A  constant  current  density  was  assumed.  Figure  8  identifies  the  t'oni 
regions  of  interest.  Regions  I,  III,  and  IV  have  a  permeability  equal  to  the 
permeability  of  vacuum,  Po,  which  for  convenience  will  be  labeled  p^.  The  problem' 
spherical  symmetry  suggested  that  a  spherical  coordinate  system  such  as  that  shown 
in  Figure  1  be  used  in  the  problem  solution. 

The  form  of  ^he  potential  in  each  of  the  regions  (1  through  IV)  was  determine^ 
from  the  solution  of  the  vector  Laplace's  equation  in  each  region.  These  magneto¬ 
static  vector  potentials  in  regions  I  through  IV  are: 
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Figure  8  -  Ferromagnetic  Spherical  Shell  Surrounded  by  an 
Infinitesimally  Thin  Current  Band 
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where  for  the  A^j  component  B  ^  =  0,  because  at  r  =  0  the  potential  must  be  finite, 
and  for  the  A^IV  component  A  ^  =  0,  because  as  r  approaches  infinity  the  potential 
must  remain  finite. 

At  each  interface,  the  basic  laws  of  magnetostatics  in  Equations  (3a)  and  (?M 
reduce  to  boundary  conditions  on  B  and  H  that  can  be  used  to  evaluate  the  six  con¬ 
stants  in  Equations  (59).  From  Equation  (18a),  the  normal  component  of  B  across  each 
boundary  must  be  continuous,  i.e.,  (B^-B^)  •  n ^  =  0  where  the  quantity  "n^7  is  the 
unit  outward  normal  to  the  surface.  This  provides  the  following  boundary  conditions 
which  must  be  satisfied  by  the  solution  in  Equations  (59)  for  each  region. 
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The  normal  component  of  the  magnetic  field  is  expressed  in  terms 
vector  potential  as 


B  =  (VxA) 
r  r 


B  =  - t — q  'Tq  (sin  9  A.) 

r  r  sin  0  30  \p 


where  B  =  V  x  A  = 


2  •  n 

r  sm  0 


3_ 

3r 


e9r 


3_ 

30 


0  0 


e  ,  r  sin  0 


3_ 


A,r  sin  0 

♦ 


However,  because  the  vector  potentials  in  each  region  are  functions  of  P‘ 


can  simplify  Equations  (60)  to  constraints  on  A 


A  =  AII  at  r  = 


AII  =  AIII  at  r  =  R2 


AIII  AIV  at  r  R3 


The  second  set  of  boundary  conditions  is  obtained  from  Equation  (18 
tangential  component  of  H  across  each  boundary  must  satisfy  the  r  ’at  ion. 
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where  (which  equals  J(0))  is  the  real  surface  current  density  in  the  limit  of 
vanishing  width  between  the  two  regions.  Using  the  relationship  B  =  pH,  Equation 
(63)  may  be  expressed  as 


B02  B01 


=  J(0) 


Referring  to  the  curl  in  Equation  (61),  we  can  write  B„  as 


B0  =  (VxA)q=  “r  h  [r  AJ 
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The  general  expressions  for  the  potentials  in  each  region  (Equations  (59))  are 

then  substituted  into  the  boundary  conditions  (Equations  (62)  and  (66))  and  solved 

for  the  constants  Api  and  B  ^  There  are  six  algebraic  equations  with  six  unknowns 

and  the  potential  in  each  region. can  then  be  specifically  determined.  The  six 

boundary  value  equations  that  must  be  solved  for  the  coefficients  are  given  next 

(where  the  index  p  is  odd  only  and  understood  to  take  on  values  from  1  to  °°) .  It  is 

noted  that  the  current  J^(0)  must  be  expanded  into  a  set  of  associated  Legendre 

functions  in  order  to  evaluate  the  constants  A  .  and  B  ..  The  detailed  expansion  is 

pi  pi 

derived  in  the  section  of  Reference  7  entitled  "Expansion  of  the  Current  (J  (0))  in 
Associated  Legendre  polynomials." 
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The  solution  of  these  equations  to  obtain  B  ^  in  terms  of  known  quantities  is  per 
formed  in  Appendix  A  of  Reference  7.  In  summary: 
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2  '  “2  / 


(68d) 


R  P+R 


P  (0)  = 


V„ (0) 
_  1  P 


Pp(cos  0)R3(p-1)(2p+l) 


(68e) 


45 


The  numerical  values  for  the  other  five  coefficients  can  be  obtained  from  the 
following  equations: 


V  M  +  Jp'  <e> 


(69a) 


A  -  =  j'  (6) 

p3  p 


Ap2  ‘  M  Bp2 


Apl  "  Ap2  +  Bp2Rl 


-(2p+l) 


(69b) 


(69c) 


(69d) 


!p*  ‘  Ap3R3<2P+1>  +  Bp3 


(69e) 


Because  the  coefficients  A  .  and  B  .  can  be  determined  from  Equations  (68)  and 

pi  pi 

(69),  Equations  (59)  can  now  be  used  to  completely  specify  the  potentials  A^,  A^^, 
Ajj. j ,  and  AIV  in  regions  I  through  IV.  The  normal  (B^)  and  tangential  (BQ) 
components  of  the  magnetic  induction  in  regions  I  through  IV  can  be  determined  by 
using  Equations  (61)  and  (65)  as: 


oo 

B0!  =  "  ^  (P+1)  (Apir(P_1))  Pp  (GOS  9) 


(70a) 


B 


rl 


1 

sin  0 


£  [APir(p_1)]  Je  (sine  (cos  e)) 


(70b) 
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Oil 


E  [Ap2(p+l)r(P-1>  -  PB^r-(P+2)l  P*  (cos  6) 


p2  J  P 


(70c) 


•p=L 


rll 


£,  iib  [vt<P"1>  +  V~<P+2>]  h  (aI“ 6  plP  (c°s  6>)  OM> 

P-1 


Bill 


L  ^Ap3(P+l)r(p-1)  -  pBp3r  (P+2)J  Pp  (cos  6) 


(70e) 


Bmi  =  ^  Titr  [Ap3r(P_1)  +  BP3r_(p+2)]  le  (sin  9  pp  (cos  9))  (70f) 


B 


u  IV 


£  0>V<p+2>)  pp <co 

0=1 


s  0) 


(7  Og) 


r  IV 


(JO 

^  Ti to  [Bp4r_(p+2)]  le  (sin  9  plp  (cos  e)) 


(7  Oh) 


In  Appendix  C  of  Reference  7,  the  coefficients  Ap^  (i  =  1,2,3)  and  Bp3 
(i  =  2,3,4)  for  the  vector  potentials  for  this  ferromagnetic  shell  problem  were 
shown  to  reduce  to  the  potentials  in  the  two  regions  of  the  simple  current  band 
problem  when  the  permeability  of  the  ferromagnetic  shell  approaches  that  of  the 
surrounding  medium  This  showed  that  the  solutions  of  that  ferromagnetic  current 

problem  had  the  correct  mathematical  form. 


47 


SPHERICAL  SHELL  WITH  INTERNAL  AND  EXTERNAL  INFINITESIMALLY 
THIN  SPHERICAL  CURRENT  BANDS 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic  spherical 
shell  of  outer  radius  R^ ,  inner  radius  R£ ,  having  a  homogeneous  permeability  p9, 
surrounded  by  an  infinitesimally  thin  current  band  of  radius  R^  having  a  current 
density  and  surrounding  an  infinitesimally  thin  current  band  of  radius  having 
a  current  density  J^.  A  constant  current  density  is  assumed  for  both  bands.  Figure 
9  identifies  the  five  regions  of  interest.  Regions  I,  II,  IV,  and  V  have  a 
permeability  equal  to  the  permeability  of  vacuum,  pQ,  which  for  convenience  will  be 
labeled  p^.  The  problem's  spherical  symmetry  suggests  that  a  spherical  coordinate 
system  such  as  that  shown  in  Figure  1  be  used  in  the  problem  solution.  The  govern¬ 
ing  differential  equation  for  A  when  homogeneous  and  linear  materials  are  considered 
is,  from  Equation  (17), 


4*  A  =  -pj 


(71) 


From  Equation  (22),  we  see  that  the  elemental  vector  potential  ctA  due  to  a 
current  element  Jdv  is  in  the  same  direction  as  J.  It  is  well  known  from  this  that 
the  lines  of  the  magnetic  vector  potential  A  are  circles  centered  about  the  coil  or 
loop  axis.  The  magnitude  of  A  along  such  a  circle  is  constant,  which  means  that  A 
is  a  function  of  the  spherical  coordinates  r  and  0  only.  Therefore,  we  know  in 
advance  for  this  problem  that  A^  is  the  only  component  of  A  existing  at  the  field 
point.  The  infinitesimally  thin  bands  of  current  shown  in  Figure  9  have  only  an 
azimuthal  or  ip  component,  which  is  a  function  of  r  and  0,  and  lie  on  the  boundaries 
between  regions  I  and  II  (i.e.,  r  =  R^)  and  between  regions  IV  and  V  (i.e.,  r  =  R^). 
These  currents  can  be  expressed  as: 


J1  =' 


vV0)’ 


if  0  <  9^  or  0  >  0£ 
if  0:£  0  _<  ©2 


(72) 
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0 


if  0  <  0^  or  0  >  0^ 


(73) 


i 

ij. 


k 


► 


I 

i 


I 

k 


V*<e) 


if  e'  <  e  <  e'' 


Therefore,  Equation  (71)  has  only  an  azimuthal  component  and  can  be  expressed  as: 


(r,6)  =  0  (in  regions  I  through  V) 


(74) 


When  the  vector  Laplacian  is  expanded  in  spherical  coordinates.  Equation  (74) 

can  be  written  as 


_2  1_  A^P  cot  9 

t  3r  t2592  r2  38 


2 


r 


0 


0  in  regions  I  (75) 
through  IV 


To  solve  Equation  (75),  we  follow  the  procedure  given  on  page  17.  The  general  solu¬ 
tion  of  Equation  (75)  in  regions  I  through  IV  may  be  formed  from  the  product  of  the 
solutions  in  Equation  (31)  which  yields 


A  =  R(r)0  (0) 


R  (r)0  (6) 
P  P 


(76) 


oe 

L 

p=i 


(c  p1 
V  p  p 


(cos  0) 


+  D  o1 

P  P 
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In  the  spherical  case,  associated  Legendre  functions  of  the  second  kind  are 
infinite  at  cos  0  =  ±  1,  and  thus  cannot  be  included  when  the  region  under  considera¬ 
tion  incLudes  the  symmetry  axis.  Therefore,  the  constant  D  must  be  set  equal  to 
zero,  and  Equation  (76)  reduces  to 


(cos  0) 

P 


(77) 


where  A  =  a'c  ,  and  B  =  b' C  . 

P  P  P  P  P  P 

The  form  of  the  potential  in  each  of  the  regions  (I  through  V)  is  determined 
from  Equation  (77).  These  magnetostatic  vector  potentials  in  regions  I  through  V 
are : 


Ar  =  A. 

I  ijj 


T  =  I]  A  -,rP  P1  (cos 

1  Lpl  J  ? 


(78a) 


ali  ■  \n  ■  £  V"  +  PP  <cos  6) 

p=l  L  r  J 


where,  for  the  A^  component,  B  ^  =  0,  because  at  r  =  0  the  potential  must  be  finite 
and,  for  the  A^v  component,  A^,.  =*  0,  because  as  r  approaches  infinity  the  potential 
must  remain  finite. 

At  each  interface,  the  basic  laws  of  magnetostatics  in  Equation  (2a)  reduce  to 
boundary  conditions  on  B  and  H  that  can  be  used  to  evaluate  the  eight  constants  in 
Equations  (78).  From  Equation  (18a),  the  normal  component  of  B  across  each  boundary 
must  be  continuous,  i.e.,  (B^  -  B^)  *  n^  =  0  where  the  quantity  n^  is  the  unit  out¬ 
ward  normal  to  the  surface.  This  provides  the  following  boundary  conditions  which 
must  be  satisfied  by  the  solution  in  Equations  (78)  for  each  region. 


Brl  =  BrII  at  r  "  R1 


BrII  ~  BrIII  at  r  ~  R2 


B  TTT  =  B  T,r  at  r  =  R~ 
rill  rIV  3 


(7  9) 


B  =  B  „  at  r  =  R. 
rIV  rV  4 


The  normal  component  of  the  magnetic  field  B^ 
vector  potential  as 


is  expressed  in  terms  of  the 


B  =  ( VxA) 
r  r 


„  _  1  3  (sin  6  A,.) 

Br  “  7  sirHe  36  ^ 


(80) 
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where  B  =  V  x  A  = 


2  •  Q 

r  sm  0 


en  r  e  .  r  s  in 

0  ip 


3  _ 

30  dip 


A.r  sin  0 


However,  because  the  vector  potentials  in  each  region  are  functions  of  P  (cos  0) 
we  can  simplify  Equation  (79)  to  constraints  on  A,,, 


Aj  =  Ajj  at  r  =  R 


(81a) 


AII  AIII  at  r  R2 


AIII  AIV  at  r  =  R, 


V  =  \  at  r  = 


The  second  set  of  boundary  conditions  is  obtained  from  Equation  (18b).  The 
tangential  component  of  ll  across  each  boundary  must  satisfy  the  relationship 


nl2  X  (H2  Hl)  Js 


where  .1  (which  equals  J(n))  is  the  real  surface  current  density  in  the  limit  of 
vanishing  width  between  the  two  regions.  Using  the  relationship  B  =  i/li,  Equation 
(82)  may  be  expressed  as 


=  JO) 


(83) 


B62  hi 
y2  “  Mi' 


Referring  to  the  curl  in  Equation  (80),  we  can  write  B„ 

t) 


as 


(84) 


From  Equations  (82)  through  (84)  the  tangential  components  in  regions  I 
through  V  must  satisfy  the  relationships: 


Lil 

r  3r 


(rAII*  +  u1  r  3r  (rAi} 


=  J^(0)  at  r  = 


(85a) 


1_  1  3__ 
)J  2  r  3r 


(rAIII> 


7  h  (CAXI> 


=  0  at  r  =  R, 


(85b) 


i_  I 

r  3r 


(rAIV)  + 


1  JL 

r  Sr 


(rAm) 


0  at  r  =  R^ 


(85c) 


1_  1  3_ 
yj_  r  3r 


(r\) 


13  f  K  \ 

7  37  (rAiv) 


J2(0)  at  r  = 


(85d) 


The  general  expressions  for  the  potentials  in  each  region  (Equations  (78))  are 
then  substituted  into  the  boundary  conditions  (Equations  (81)  and  (85))  and  solved 
for  the  constants  ^  and  B  There  are  eight  algebraic  equations  with  eight 
unknowns,  and  the  potential  in  each  region  can  then  be  specifically  determined.  The 
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eight  boundary  value  equations  that  must  be  solved  for  the  coefficients  are  given 
next  (where  the  index  p  is  odd  only  and  understood  to  take  on  values  from  1  to  ®). 
It  is  noted  that  the  current  J^(0)  must  be  expanded  into  a  set  of  associated 
Legendre  functions  in  order  to  evaluate  the  constants  A  ±  and  B  i<  The  detailed 
expansion  is  derived  in  the  section  of  Reference  7  entitled  "Expansion  of  the 
Current  (J^(0))  in  Associated  Legendre  Polynominals . " 

4pX  ’  [V?  +  Bp2Rl<P+1>]  <8Ra> 


v;<p+1>]  ■ 


p3  2 


+  B 


p3 


-(^)j 


(86b) 
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B  ,R  (p+1) 
p4  3  ... 


(86c) 
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(8  6d) 
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1LAP2(P+1)R1 


(p-D 
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(p+1)  A  R(p_1) 
pi  1 


Jpl(Q> 

p1 (cos e  ) 
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(86e) 


-  ^[Ap3<P+1>R2P_1>-PBp3RItP+2;‘J  +  yVP+1)R2','1)-PBp2R;l[>+2>]  ’  0  <««> 

"  S7  [V<P+1)R3<P  1>"pBp4,3  <P+2>]+  Uj[Ap3<P+1)R3P"1>-PBp3R(<P+2)]  -  0  <86s) 
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The  solution  of  these  equations  is  performed  in  Appendix  D.  In  summary: 


A 


yP2<9> 


p4  {2p+l)R^p-1)  Pp  (cos  6  ) 


5  [X] 


Vpl(6> 


p2  (2p+l)R"(p+2)PP  (cos  6) 


[Y] 


P3 


m  ->) 


p2 


p2 


1  + 


(2p+l) 


aP3  =  Bp3M  +  M[s]=  [z] 


P3 


RP[X][S.]  -  RP[X][S]rT]  -  [Y][A]RP  -  [Y]R~(p+1) 
RP  fW][T]  -  R~(p+1)  -  RP  [W] 
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(87a) 


(87b) 


(87c) 


(87d) 


(87  e) 


Because  the  coefficients  A  and  B  .  can  be  determined  from  Equations  (86)  and 

pi  pi 

(87),  Equations  (78)  can  now  be  used  to  completely  specify  the  potential  s  AJf  An, 
Am,  and  A^,  in  regions  I  through  IV.  Then  the  normal  (B^)  and  tangential  (Bq) 
components  of  the  magnetic  induction  in  regions  I  through  IV  can  be  determined  by 
using  Equations  (80)  and  (84),  to  be: 


E  [(p+DAplr 


(p-15Jp^(cos  0) 


(88a) 


.(p-1) 


rl  sin  6  \  PI 

p^l 


£  (AP!r(P_1))  ^  [Sin  6  Pp  (C°S  9)] 


(88b) 


;I  -  Z  [-<p+u  +  pSp2'‘(P+2)]  0> 


(88c) 


rll  sin 
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(88d) 
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••III  ,  •"  P5  J  P 


(88e) 
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P  (cos  0) 
P 
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In  Appendix  D,  the  coefficients  A  ^  (i  =  1,2, 3, 4)  and  B^  (i  =  2, 3, 4, 3)  for 
the  vector  potentials  for  the  present  ferromagnetic  shell  problem  reduce  to  the 
potentials  in  the  two  regions  of  the  simple  current  band  problem  when  the  permeabil¬ 
ity  of  the  ferromagnetic  shell  approaches  that  of  the  surrounding  medium  Uj. 

This  shows  that  the  solutions  of  the  above  ferromagnetic  current  problem  have  the 
correct  mathematical  form. 

SOLUTIONS  FOR  PROLATE  SPHEROIDAL  BODIES 

SOLID  PROLATE  SPHEROID  OR  PROLATE  SPHEROIDAL  SHELL  IN  A  UNIFORM 
FIELD  OF  ARBITRARY  DIRECTION 

Important  problems  relating  to  determining  the  magnetic  induction  for 
spheroidal  ferromagnetic  bodies  in  an  external  magnetic  field  have  not  been  widely 
reported  on  in  the  literature  or  in  text  books.  The  solutions  for  these  types  of 
boundary  value  problems  can  be  obtained  by  using  a  procedure  similar  to  that  used 
for  spherical  bodies"*’  ’  by  determining  the  solution  to  Laplace's  equation  for  the 


59 


magnetic  scalar  potential.  Constant  external  field  problems,  solid  and  shell,  were 
solved  and  programmed  on  the  digital  computer  by  Nixon  of  the  Center,  for  the  case 
of  an  arbitrarily  oriented  external  magnetic  field. ^  The  solutions  found  in 
Reference  6  were  presented  in  Cartesian  coordinates.  The  problem  of,  for  instance, 
finding  the  magnetic  induction  for  an  infinitesimally  thin  current  band  surrounding 

Q 

a  spheroidal  shell  can  be  generalized  to  include  an  external  magnetic  field. 

Linear  superposition  may  be  applied  to  find  the  solution  in  this  case.  Therefore, 
in  Appendix  E  the  Cartesian  expressions  were  converted  to  prolate  spheroidal 
coordinates  to  be  compatible  with  other  problem  solutions  in  this  section  of  the 
report . 

SOLID  PROLATE  SPHEROID  SURROUNDED  BY  AN  INFINITESIMALLY 
THIN  SPHEROIDAL  CURRENT  BAND 

8  9 

For  the  case  of  the  prolate  spheroidal  bodies  *  the  equations  given  in  the 
Basic  Equations  section  of  this  text  apply.  The  governing  differential  equation 
for  A  when  homogenous  and  linear  materials  are  considered  is  from  Equation  (17). 


2$  A  =  -pJ 


(89) 


where  the  general  expression  in  prolate  spheroidal  coordinates  for  a  current  density 
is 


J 


+  J 


0 


+  J 


'P  IP 


(90) 


As  previously  discussed  (see  page  16),  because  the  current  has  only  an  azimuthal  or 
ip  component,  A  has  only  an  A^  component.  For  the  spheroidal  problems  considered  in 
this  report,  the  current  band  is  assumed  to  be  infinitesimally  thin  and  the 
governing  differential  equation  for  A  in  each  region  can  be  expressed  as 
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1  9(sinh  n  A^) 

.3 

1 

d(sin  0  A.) 

'P 

sinh  n  dp 

+  JQ 

sin  0 

90 

(91) 


Using  the  method  of  separation  of  variables,  the  solution  to  Equation  (91)  is 


T.  [a  P'*‘(cosh  n)  +  B  Q^(cosh  rm 

P=i  L  p  p  J 


X 


p\cos  0)  +  B/ 
P 


Q1(cos  0) 
P 


(92) 


where  Pm  and  Qm  are  the  associated  Legendre  functions  of  the  first  and  second  kind, 

P  P 

respect ively . 

For  the  prolate  spheroidal  system,  the  associated  Legendre  functions  of  the 
second  kind  are  infinite  at  cos  0  =  ±1,  and  as  such  cannot  be  included  in  a  general 
solution  for  a  given  region  which  includes  0=0,  or  0  =  Tt.  Therefore,  in  our  case, 
the  constant  ?>'  is  set  equal  to  zero.  Equation  (92)  reduces  to 


A .  = 


£  l  i  p 

p=i 


k.,pl(cosh  n)  +  k^Q^(cosh  ri) 


2  P 


J  P^  (cos  0) 


(93) 


where  k^  and  k ^  are  constants  (k^  =  AA  ,  k^  =  BA^  ).  When  the  substitutions  £ 
cosh  r|  and  v  =  cos  0  are  made  in  Equation  (93),  A,  can  be  expressed  as 


OO 

\  =  £  [v J(C)  +  k2Qp(o]  pj 


(V) 


(94) 
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This  is  the  general  form  of  the  psi  (i^)  component  of  the  vector  potential  that  will 
be  used  to  determine  the  potentials  A^  in  each  region. 

The  problem  of  a  solid  ferromagnetic  prolate  spheroid  surrounded  by  an 
infinitesimally  thin  prolate  spheroidal  current  band  shown  in  Figure  10  was  solved 
by  Purczynski.^ 

In  this  case,  the  permeability  of  the  solid  spheroid  is  and  the  boundary  of 
the  body  is  determined  by  n  =  n^  =  constant.  The  permeability  yQ  of  a  vacuum  that 
is  external  to  the  spheroid  is  denoted  by  y  The  current  band  which  lies  in  the 

boundary  between  regions  II  and  III  is  denoted  by  p  =  p^  =  constant,  and  the  con¬ 
stant  current  density  flowing  in  the  band  is  J. 

For  completeness,  Purczynski's  work^  is  presented  in  this  text  in  our  nota¬ 
tion.  The  form  of  the  components  of  the  vector  potential  in  regions  I  through 
III  is  determined  from  Equation  (94).  These  magnetostatic  vector  potentials  in 
regions  I,  II,  and  III  are: 


V  •  £  [AppP1({)]  pP1(v>  (95a> 

oo 

v  ■  E  [Vp«>  +  Vp(w]  PJ(V)  (95b> 

P=1 

Vu-  £  pp(u)  (95c> 

p=i 

where  £  =  cosh  p  and  v  =  cos  0. 

We  note  two  constants  were  set  equal  to  zero  because  the  potential  must  be 
finite  in  each  of  the  regions  I  through  III  and  approach  zero  as  E,  -*■  00  in  region  III. 
The  remaining  constants  are  determined  from  the  boundary  conditions 
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(96a) 
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(96b) 


nil 


=  B 
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at  n  =  nn 


B9iii 


=  j^(e)  at  n  =  n2 


0  at  n  =  r\ 


Equations  (96)  can  be  written  at 
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The  general  expressions  for  the  potentials  in  each  region  (Equations  (95))  are  then 
substituted  into  the  boundary  conditions  (Equations  (97))  and  are  solved  for  the 
four  constants.  Because  there  are  four  equations  with  four  unknowns,  the  potential 
in  each  region  can  be  determined.  The  four  boundary  value  equations  are  presented 
below.  The  index  p  in  the  summation  sign  has  both  even  and  odd  values  and  takes  on 
values  from  1  to  °°.  It  is  noted  at  this  point  that  the  current  density  J^(0)  must 
be  expanded  into  a  s’et  of  associated  Legendre  functions  to  evaluate  the  vector 
potential.  The  detailed  expansion  is  presented  in  Reference  8. 


yj  -  [vp(y +  vP1(i;i)KM 


(98a) 


[Bppp({2)  +  Vp  («2)K(u) 


Vp  (52)]Pp<V) 


(98b) 


-ft) 


1  8 


ST 


[yj«>]  pJ(v> 


35 


^2‘1^[BpPp(E)  +  Vp*5’]  Pp 


(V) 


Jp(6) 


(98c) 
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«>  +  Vp(C> 


3 

3C 


(«2-1)'[Vp<5>]  Pp<v) 


5-5, 


(98d) 


If  we  make  the  following  substitutions 

pVj  =  [fe2-l)*  Pp(€)j  (99a) 

Qp(C)  =  [k2-if  Qp(o]  (99b) 


and  perform  simple  algebraic  manipulation,  the  four  boundary  conditions  can  be 
simplified  to: 


Vp({l)  ■  Vp(5i)  +  Vp(5i) 


(100a) 


BpPp(52)  +  Vp(?2^  "V^) 


(100b) 


J 

j 
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We  note  that  when  is  allowed,  in  the  limit,  to  approach  this  solution  reduces 
to  that  of  a  current  band  in  vacuum.  We  also  note  that  the  solutions  are  not  in  the 
identical  form  of  those  given  in  Reference  10. 
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B0III  '  ^ 


E  yr^  ^  r'l]  (vH  P>(v) 

P=1  a(c-v)  L 


( 1 02  f ) 


PROLATE  SPHEROIDAL  SHELL  SURROUNDED  BY  AN  INFINITESIMALLY  THIN 
SPHEROIDAL  CURRENT  BAND 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic  prolate 
spheroidal  shell  of  homogeneous  permeability  surrounded  by  an  infinitesimally 
thin  prolate  spheroidal  current  band  of  constant  current  density  J.  The  geometry 
of  the  problem  suggests  that  a  prolate  spheroidal  coordinate  system,  as  shown  in 
Figure  2,  should  be  used  in  the  solution.  Figure  11,  a  cross  section  of  the 
problem  geometry,  identifies  the  four  regions  of  interest.  The  boundaries  of  the 
prolate  spheroidal  shell  are  determined  by  I  =  =  constant  and  n  =  =  constant. 

The  constant  current  lies  in  the  boundary  n  =  =  constant.  Regions  I,  III,  and 

IV  have  a  permeability  labelled  Ampere's  law  states  that 


V  x  H  *  J 


(103) 


and,  because  7  *  B  =  0,  the  induction  B  must  be  the  curl  of  some  vector  field  A. 
The  governing  differential  equation  for  A,  when  homogeneous  and  linear  materials 
are  considered,  is,  from  Equation  (17), 


(104) 


The  general  expression  in  prolate  spheroidal  coordinates  for  a  current  density 


J  =  J  e  +  j  ef  +  J  e , 

n  n  oo  4’ 


(105) 
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In  the  problem  presented  herein,  the  current  density  has  only  a  psi  (^)  component 
fj^(0)  J ,  which  means  that  the  vector  potential  has  only  a  psi  component  A^e^. 
The  vector  potential  A  =  is  a  function  of  the  prolate  spheroidal  coordinates 


H ,  and 


i.e. 


[VVn’0)] 


The  constant  current  density,  which  lies  on  the 


boundary  between  regions  III  and  IV,  can  be  expressed  by  the  function 


(o 


J  =  < 


if  0  <  6^ or  0  >  ©2 


Ve)S 


if  0  <  e2 


(106) 


where  (0)  is  equal  to  a  constant  J  along  n  =  0-  ^or  0^  £  0  £ 

Therefore,  Equation  (104)  has  only  an  azimuthal  component  and  can  be  expressed 

as 


4*  A^(n,6)  =  0 


/ in  regions  I \ 
\through  IV  / 


(107) 


When  the  psi  component  of  the  vector  Laplacian  X$X  A(  is  expanded  in  prolate 

T 

spheroidal  coordinates.  Equation  (107)  can  be  expressed  as  (see  Appendix  A) 


3 

l  ‘ 

)(sinh  n  A^)  ' 

3 

_1_  - 

r 

l 

3 (sin  0  A^)' 

_  sinh  n 

3n 

+  30 

sin  0 

30  J 

(in  regions  I  through  IV) 


(108) 


Applying  the  method  of  separation  of  variables,  let  us  assume  that  A^  can  be 
expressed  as  the  product  of  two  functions 


H(cosh  n)  "(cos  0) 


(109) 


where  H(cosh  f|)  is  a  function  of  cosh  ri  only  and  G(cos  0)  is  a  function  of  cos  0 
only.  Substituting  this  form  of  the  component  of  the  vector  potential  A  into  Equa¬ 
tion  (108)  , ,  we  have,  after  separation  of  variables, 


} 

i +  coth  n  ^  -  (p(p+1)  +  ~ .  T?  ") H 

1  \  sinh  n/ 


(110a) 


i-f  +  cot  e  jf  +  /P(p+i) 


--V-) 

sin  0  / 


G  =  0 


(110b) 


where  the  separation  constant  is  p(p+l)  and  p  is  an  integer  from  one  to  infinity. 
It  is  well  known  that  differential  equations  of  the  form 


d‘H  ,  dH 

— -  +  coth  n  - — 

dp2  dri 


have  the  general  solution  of  the  form 


=  C,Pm(cosh  n)  +  C0Q"'(cosh  f|) 

1  p  2  p 


p(p+l)  +  B  2  H  =  0 
s in  h  n  / 


(111a) 


(111b) 


where  and  C ^  are  constants.  It  is  known  that  a  differential  equation  of  the 
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hV 

— +  cot0 
<39 


p(p+l) 


2 

m 

sin^  0. 


/ 

G 


0 


(112a) 


has  the  general  solution  of  the  form 


G  =  C-,Pm(cos  0)  +  C,Qm(cos  0) 
J  p  4  p 


(112b) 


where  and  are  constants.  The  associated  Legendre  functions  pm  and  are  of 

the  first  and  second  kind,  respectively.  Comparison  of  Equations  (110),  (111),  and 

2 

(112)  shows  that  in  Equations  (111)  and  (112),  m  =1.  This  requires  that  m  always 
equals  unity.  The  solutions  of  Equations  (110a)  and  (110b)  are  expressed  as 


H(cosh  n)  =  A  P1(cosh  n)  +  B  Q1(cosh  n) 
P  P 


(113a) 


G(cos  0) 


A 


P^(cos  0)  + 
p 


(cos  0) 


(113b) 


The  general  solution  of  Equation  (108)  may  be  formed  from  the  product  of  solutions 
to  Equations  (113a)  and  (113b)  which  yield 


A  =  H(cosh  n)  G(cos  0)  =  H  (cosh  n)  G  (cos  0) 

V  P  P 


(114) 


oo  ^ 

=  E  |A  P1(cosh  n)  +  B  Q1(cosh  r))| 
p=l  L  p  P  J 

(115) 


x  ^A'  Pp(cos  9)  +  Br  Q^(cos  0)J 


For  the  prolate  spheroidal  system,  the  associated  Legendre  functions  of  the  second 
kind  are  infinite  at  cos  0  =  ±1  and,  as  such,  cannot  be  included  in  a  general  solu¬ 
tion  for  a  given  region  which  includes  0  =  0  or  0  =  tt.  Therefore,  in  our  case,  the 
constant  is  set  equal  to  zero.  Equation  (115)  reduces  to 


OO 

A.  =  X  [~K  P^cosh  n)  +  K.Q1(cosh  n)l  P^cos  0) 
V  .  Lip  2  p  J  P 


(116) 


where  and  are  constants  (K^=AAX  ,  K2=AB/ ) .  When  the  substitutions  £  «  cosh  r 
and  v  =  cos  0  are  made  in  Equation  (116),  A,  can  be  expressed  as 


OO 

\  •  ?!  [K1PP(S)  +  K2QP<f;)]  Pp' 


(117) 


This  is  the  general  form  of  the  psi  component  of  the  vector  potential  that  will  be 
used  to  determine  the  potentials,  A  ,  in  each  region. 

The  form  of  the  component  of  the  vector  potential  A^  in  regions  I  through  IV  is 
determined  from  Equation  (117).  These  magnetostatic  vector  potentials  in  regions 
I  through  IV  are: 
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v  ■  P4  [vH  pJ  (“> 

-  £  [Vp(5)  +  V^0]  ?1P<V> 

(118) 

OO 

a)IjTTT  =  £  I’d  p*(0  +  E  oJ(5)l  pj-(v) 

V111  P=1  L  p  P  P  P  J  p 

OO 

Vv  ■  £  [Vp<5)]  PP<',) 


Bi'i  nuso  the  potential  must  be  finite  in  each  of  the  regions  I,  II,  and  III  and 
approach  zero  as  (  +  in  region  IV,  the  following  constants  were  set  equal  to  zero. 


For  A  ,  the  constant  associated  with  Q^(£)  P^(v)  was  set  equal  to  zero 

'•P 1  P  P 


hoi  au  se 


Qp(£)  -*•  ®  at  £  =  +1 


J.  For  A  ...  the  constant  associated  with  P^(£)  P^(v)  was  set  equal  to  zero 
j  viv  p  P 

because  p  (■*)  ■*  ®  as  £  ■*  °° . 

P 

(We  note  0^(£)  -*■  0  as  E,  -*  <*>) 

P 


The  constants  A  ,  B  ,  C  ,  D  ,  E  ,  and  F  are  to  be  determined  from  the  boundary  con- 
P  P  P  P  P  P 

ditions.  At  each  interface,  the  basic  laws  of  magnetostatics  (Equations  (3a)) 
reduce  to  boundary  conditions  on  B  and  H  that  can  be  used  to  evaluate  these  six 
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constants. 


The  normal  component  of  B  across  each  boundary  must  be  continuous,  i.e., 
•  n12  =  0  where  the  quantity  is  the  unit  outward  normal  to  the  surface. 
This  provides  the  following  boundary  conditions  which  must  be  satisfied  by  the  solu¬ 
tions  given  in  Equation  (118)  for  each  region. 


Bni  =  Bnn  at  n  =  ni 


(119a) 


Bnn  -  B  XI!  at  n  =  n2 
n 


(119b) 


Bnni  ‘  Bniv  *t  "  '  "3 


(119c) 


The  eta  (r|)  or  normal  component  of  the  magnetic  field  (B^)  is  expressed  in  terms  of 
the  vector  potential  as 


(120) 


where  B  =  V  x  A 


_ 1 _ 

a(sinh2  n  +  sin2  o) (sinh  n  sin  0) 


x 


(Note:  Above  equation  continued  on  next  page). 


e'  (sinh^  n 
n ' 


+  sin  0 


)  «»( 


Isinn  n  +  sin  Q) 


sinh  n  sin  6 


s_  s_  a_ 

3n  30  3ip 


0 


0 


sinh  r)  sin  0 


and 


£,  =  cosh  n, 


{  2  2 
isinh  n+  sin  0 


V  =  cos  0,  *  a  sinh  ri  sin  0 

However,  because  the  vector  potentials  in  each  region  are  functions  of  P^(v), 
can  simplify  Equation  (119)  to  constraints  on  at  the  interfaces: 


\l  =  \ll 


\ll  ~  AipIII 


\lll  ~  \lV 


at  n  =  n1 


at  n  *  n2 


at  n  =  n3 


The  second  set  of  boundary  conditions  states  that  the  theta  (0)  or  tangential 
ponent  of  H  across  each  boundary  must  satisfy  the  relationship 


we 


(121a) 


(121b) 


(121c) 

com- 


(122) 
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where  Jg  (which  equals  J^(0)e^  is  the  real  surface  current  density  in  the  limit  of 
vanishing  width  between  the  two  regions.  Using  the  relationship  B  =  pH,  Equation 
(122)  can  be  expressed  as 


V0> 


(123) 


Referring  to  the  curl  in  Equation  (120),  we  can  write  B„  in  the  form 


Ba  =  (VxVe 


-l 


ele3 


(124) 


From  Equations  (123)  and  (124)  the  tangential  components  of  B  in  regions  I  tnrough 
IV  must  satisfy  the  relationships: 


(123a) 


(125b) 


=  Z  JpO) 


oo 

£ 

P=1 


KG  P1(v) 
P  P 


(125c) 


The  general  expressions  for  the  potentials  in  each  region  (Equation  (118))  are  then 
substituted  into  the  boundary  conditions  (Equations  (121)  and  (125))  and  solved  for 
the  six  constants  (A  ,  B  ,  C  ,  D  ,  E  ,  and  F  ).  Because  there  are  six  equations 

p  p  p  p’  p’  p  n 

with  six  unknowns,  the  potential  in  each  region  can  be  determined.  The  six  boundary 
value  equations  are  presented  below.  The  index  p  in  the  summation  sign  has  both 
even  and  odd  values  and  takes  on  values  from  1  to  °°.  It  is  noted  at  this  point  that 
the  current  density  .1^(0)  must  be  expanded  into  a  set  of  associated  Legendre 
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functions  to  evaluate  the  constants  in  the  vector  potential  (Equation  (118)) 
detailed  expansion  is  presented  in  Appendix  B  of  Reference  8.  The  six 
expressions  for  the  boundary  conditions  are: 


.  The 


(126a) 


(126b) 


(126c) 


(126d) 
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PL(V)  =  F 
P  ■ 


PQJ(E3)'p(V) 


If  we  make  the  following  substitution 


Pp<^> 


d 

dC 


(^) 


P1(Q 

P 


Qp(^> 


d 

d£ 


(52-l) 


and  perform  simple  algebraic  manipulations,  the  six  boundary  conditions  can 
plified  to: 


Vp(Fl)  ■  Vp(-l)  +  CpQp(h) 


(126e) 


(126f) 


(127) 


(128) 

be  sim- 
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The  solution  of  these  six  simultaneous  equations  to  obtain 


quantities  gives: 


in  terms  of 


kWJp  Pp<H 

K) 

-p(y 

p  k' 

(wwpp(H 

)^( 

w<h; 

2) 

M- 


i(h) 

Wpp({2)  +  <(«2) 


known 


(130a) 


(130b) 


(130c) 


(130d) 


(130e) 
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IP1 

TI  _  'ppV7 


WPp(?2)  +  ^2) 


The  numerical  values  for  the  other  five  coefficients  can  be  obtained  from  the 
following  equations: 


C  =  J  +  E  [2] 
P  P  P 


B  =  TxlC 
P  P 


D  -j1 
P  P 


*  +c 

p  p  P  Pi(Si) 


Vpfe) 

<  N 


Because  the  six  coefficients  can  be  determined  for  a  specified  problem  from  Equa 
tions  (13C)  and  (131),  the  potentials  A^T,  a^ju»  an<^  Aipiv*  re6i°ns  1 

through  IV  can  be  completely  determined.  The  normal  (B  )  and  tangential  (BQ)  to 


the  surface  n  =  constant  (or  5  =  constant)  components  of  tne  magnetic  induction  in 
each  region  I  through  IV  can  be  determined  by  using  Equations  (120)  and  (124),  to  be: 


ol 


=  -  E 


^  a(^-v2) 


1  8 
h  35 


fe2-l)  A  pi(OpJ(v) 


P  P  P 


ni 


-  E 


P=1  afe2-v2) 


3  fv.  2)\  P1(5)pJ'(v)| 

%  "3v  LU*V  '  p  p  p  J 
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=  -  E 
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!  t  2  2V* 
P=1  a(52-v2) 


+  cpoJ(e))pJ(v) 


mi 


^[(l-^^pjco  +  Cpq‘(5))pJ(v)] 


0111 


^  4i2-*2)h 


k  (»: 


pX(5)  +  E  <£(5) 
p  p  p  p 


’)pM 


Vu  E 


p-1  4c2-v2)' 


|_[(l-v2f(DppJ(«  +  EpQJ,C))  pj(v) 


(132a) 


(132b) 


(132c) 


(132d) 


(132e) 


(132f ) 
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(132g) 
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PROLATE  SPHEROIDAL  SHELL  SURROUNDING  AN  INFINITESIMALLY  THIN 
SPHEROIDAL  CURRENT  BAND 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic 
spheroidal  shell  of  homogeneous  permeability  surrounding  an  infinitesimally 
thin  prolate  spheroidal  current  band  having  a  constant  current  density  J.  Figure 
12  shows  the  cross  section  of  the  problem  geometry.  The  coordinate  system  shown 
previously  in  Figure  2  will  be  used  in  the  solution.  The  boundaries  of  the  prolate 
spheroidal  shell  are  determined  by  n  =  and  n  =  U2 .  The  steady  state  current 
lies  in  the  boundary  n  =  and  between  0^  <  0  <  As  in  the  previous  problem, 

the  constant  current  density  has  only  a  psi  component  J^(0)e^,  and  thus  the  vector 
potential  has  onl’  psi  component  A^e^ .  The  vector  potential  is  a  function  of  the 
prolate  spheroj.  coordinates  n  and  0.  The  constant  current  density  is  expressed 
by  Equation  (106)  when  the  boundary  n  is  changed  to  n  =  0^. 

The  governing  partial  differential  equation  has  only  a  psi  component  and  is 
given  by 


A  =  $5A^(r|,H)  =  0  (in  regions  I  through  IV)  (133) 


When  the  vector  Laplacian  $  A  is  expanded  in  prolate  spheroidal  coordinates. 
Equation  (133)  can  be  expressed  as  (see  Appendix  A  of  Reference  8) 


9 

r  i 

9r| 

sinh  ri 

9(sinh  n  A^) 


-h  M; 


sin  0 


9_ 

90 


(sin  u 


A,) 

J 


(134) 
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Figure  12  -  Infinitesimally  Thin  Current  Band  Surrounded  by 
Ferromagnetic  Spheroidal  Shell 


Adopting  the  following  notation 


E,  =  cosh  n,  —  cos  0  (135) 

and  following  the  logic  presented  earlier,  the  solutions  for  in  regions  I 
through  IV  have  the  general  form 

A  =  i  +  m[(o]p!(v)  (136) 

»P  P=1  L  p  p  J 

The  form  of  the  components  of  the  vector  potential  A^  in  each  of  the  regions 
I  through  IV  is  determined  from  Equation  (136) .  These  components  of  the  vector 
potential  in  each  region  are: 


V-  i  [HppP(0lpP(v) 
p=i  *-  -* 

Vt  ■  £  [vJ(E)  +  V^<£)lPp(v) 

P=1  L  -I 

Vn  ■  £  [vp<S)  +  Vp<5)]  pp(v) 
Vv  ■  £  [»P<<5>]pp<v> 


(137a) 


(137b) 


(137c) 


(137d) 


The  functions  are  the  associated  Legendre  functions  of  the  first  kind  of  degree 
1  and  order  p,  and  the  Q1  functions  are  associated  Legendre  functions  of  the 


second  kind. 


At  each  interface,  the  basic  laws  of  magnetostatics  reduce  to  boundary  condi¬ 
tions  on  B  and  H  (see  Equations  (119)  and  (122))  that  can  be  used  to  determine 
related  boundary  conditions  on  A: 


Aif>I  ~  Ai|iII  n  =  ni 


(138a) 


AipII  ~  AipIII  n  =  n2 


(138b) 


Aif>III  ~  AipIV  n  =  n3 


(138c) 
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f= 
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These  boundary  conditions  are  then  used  to  evaluate  the  constants  in  Equation  (137). 
Using  Equations  (137)  and  (138)  to  solve  for  the  coefficients  (where  the  index  p 
takes  on  all  values  from  1  to  00 )  we  get: 


h  p1(c  )px(v)  =  fi  p];(5  )  +  k  qV  )"|pJ 
ppip  Lpp1  pp1Jf 


(V) 


(139a) 
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.2  2 


L  P1(Q 
P  P 


MQ(C) 

P  P 


PX(v) 

P 


If  we  make  the 


following  substitutions 


Pp<5>  ’ 


’  If 


(t2-0Vco 


and  perform  simple  algebraic  manipulations,  the  six  boundary  conditions 


H  P1^)  =  I  P1(C1)  +  K  Q*(L) 

ppl  p  p  X  p  p  J. 


+  K  Q^C,) 
P  P  1 


ft)' 


■=—  j  H  PQ(^)  = 

U_  I  p  p  T. 


JP<e)»(f2-"2) 

PX(v) 


(139f) 


(14  0) 


(141) 


reduce  to : 


(142a) 


(142b) 
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(142c) 


I  ?la.)  +  K  qJ;(S,)  =  l/cL)  +  M  qJ(5  ) 

p  p  2  p  p  2  p  p  2  p  p  2 


(~)  U  Pfi(£  )  +  M  q4(C)  PA(S,)  +  K  Q4(5,)\ 

W]  lp p  2  P  P  2  )  \Ml)\  P  P  2  P  P  2  / 


( 14  2d ) 


L  P1(5,)  +  M  qJ;(C-)  =  N  o\y 
P  p  3  PPJ  P  P  -> 


(142e) 


(142f ) 


It  should  be  noted  in  the  above  equations  that  the  current  density  J^(0)  was 

expanded  into  a  set  of  associated  Legendre  functions  to  evalute  the  constants  in 

the  vector  potential  components  (see  Appendix  B  of  Reference  8) . 

The  solution  of  these  six  simultaneous  Equations  (142a)  through  (142f)  to 

obtain  L  in  terms  of  known  quantities  is: 

P 


L  = 
P 


+  JXQA(g  ) 
P  P  2 


[v]p\{2>  + 


(143) 
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(144) 


(145) 


(146) 


(147) 


(148) 


The  numerical  values  for  the  other  five  coefficients  can  be  obtained  from  the 
following  equat ions : 


K  =  J1  (149) 

P  P 


Hp  =  Lp[u]  (150) 


I  = 
P 


(151) 


N 

P 


+  M 

P 


(152) 


H 

P 


I  + 
P 


VW 


(153) 


The  components  of  the  potential  in  regions  I  through  IV  can  be  determined 

because  the  coefficients  H  ,  I  ,  K  ,  L  ,  M  ,  and  N  can  be  calculated  for  a  speci- 

P  P  P  P  P  P 

fic  problem.  The  normal  (B^)  and  tangential  (By)  components  (to  the  surface 

ri  =  constant  or  £  =  constant)  of  the  magnetic  induction  in  each  region  (I  through 

IV)  can  be  determined  by  using  Equations  (120)  and  (124),  to  be: 


(154a) 


0 1 1 


=  -  2 


P=1  a(42-V2) 


(^-l) 


*  (vi(E)  +  v£(°)  pJ<v> 


B 


mi 


*  lVp<0  +  v^<0 


pj(v) 


B0in 


L  P1(Q  +  M  QX(4) 
P  P  PP 


pj<v) 


(154b) 


(154c) 


(154d) 


(154e) 
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Psl  ar-v2 


x  ( L  P2(  Q  +  H  pj;(v) 

\  P  P  P  P  P 


(154f) 


0IV  = 


-  I 

p=l 


>(52-v2) 


r  a-\^y 

h  3C 


(n  q1 
\  p  p 


JO  p  (v) 

p  /  p 


(154g) 


-  X  - r  (l-v2)  (nq1(q)  p1(v) 

p-i  a(52.v2 y  3vLv  m  p  p  ;  p  j 


(154h) 


PROLATE  SPHEROIDAL  SHELL  WITH  INTERNAL  AND  EXTERNAL  INFINITESIMALLY 
THIN  SPHEROIDAL  CURRENT  BANDS 

We  now  proceed  to  solve  the  boundary  value  problem  of  a  ferromagnetic  prolate 
spheroidal  shell  of  homogeneous  permeability  y2  with  internal  and  external, 
infinitesimally  thin,  prolate  spheroidal  current  bands  of  constant  current  density 
and  J2,  respectively.  The  geometry  of  the  problem  suggests  that  a  prolate 
spheroidal  coordinate  system  as  shown  in  Figure  2  can  be  used  in  the  problem  solu¬ 
tion.  Figure  13,  a  cross  section  of  the  problem  geometry,  identifies  the  five 
regions  of  interest.  The  boundaries  of  the  prolate  spheroidal  shell  are  determined 
hy  D  =  ^2  and  H  =  r^,  constants.  The  direct  currents  lie  in  the  boundaries  n  =  , 

and  n  =  constants.  Regions  I,  II,  IV,  and  V  have  a  permeability  equal  to  vac¬ 
uum  Pq,  which,  for  convenience,  will  be  labelled  p^.  In  the  problem  presented  here, 
the  current  densities  have  only  a  psi  (iJj)  component  ^J^(0)e^j,  which  means  that  the 
vector  potential  has  only  a  psi  component  Aihe|(i.  The  vector  potential  |a  =  A^e^J is  a 


function  of  the  prolate  spheroidal  coordinates  fU^i.e.,  A^  =  A^  (H,0)J  .  The  con¬ 
stant  current  densities,  which  lie  on  the  boundaries  between  regions  I  and  II  and 
between  regions  IV  and  V,  can  be  expressed  by  the  functions 
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J 


1 


0,  if  0  <  0^  or  0  >  8^ 

Ve)  V  lf  61  i  9  i  e2 


(155) 


where  J^(0)  is  equal  to  a  constant  along  n  =  8^  for  0^  <_  0  <_  ©2 


and 


0, 


if  0  <  0  or  6  >  02 


J2  " 


(156) 


V(0)  if  0i  1  6  i  e2 


where  V  (0)  is  equal  to  a  constant  Jj  along  f|  =  for  0,^  £  0  £  0^ •  Therefore, 
Equation  (17)  has  only  an  azimuthal  or  psi  component  and  can  be  expressed  as 


4*  A  =  4*  A  (n,0)  =  0  in  regions  I  through  V 


(157) 


When  the  vector  Laplacian  4*  is  expanded  in  prolate  spheroidal  coordinates. 
Equation  (157)  can  be  expressed  (see  Appendix  A,  Reference  8) 


9n 


8  (sinh  n  A 


sinh  H 


3n 


+ 1_  r  i_ 

80  sin  0 


8  (sin  0  A 


80 


=  0 


Adopting  the  following  notation 


(158) 


£  =  cosh  n,  v  =  cos  0 

and  following  the  logic  presented  earlier,  the  solutions  for  A^  in  regions  1 
through  IV  have  the  general  form 
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\  -  +k2^K)1  pJ(v) 

D=1  -* 


(159) 


The  form  of  the  component  of  the  vector  potential  in  regions  I  through  V 
is  determined  from  Equation  (159).  These  magnetostatic  vector  potentials  in 
regions  I  through  V  are 


>1 


00  r 

- 

P=llPP 


(O 


pX(v) 

p 


(160a) 


>  II 


E  [yj«>  +  cPQi<5)  ]  PJ<V) 

p=l  L  J 


(160b) 


VII  + 

P=1 


(v) 


(160c) 


>IV 


+  v£(5>1 

p=i  '•  J 


(V) 


(160d) 


>V 


p=i  L  J 


(160e) 


Because  the  potential  must  be  finite  in  each  of  regions  I  through  IV,  and  approach 
zero  as  +  ”  in  region  V,  the  following  constants  were  set  equal  to  zero: 

1.  For  A^j.,  the  constant  associated  with  Q^(£)  Pp(v)  was  set  e9ual  to  zero 
because 


Qho 


00  at  £  =  1  (z  axis  between  +a) 
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P^(£)  ■+  ®  as  C  +  “ 
P 


(we  note  Q1^)  +  0  as  £  •+•  00  )  • 

P 

Constants  A,B,C,D,E,F,G,  and  H  are  to  be  determined  from  the 
P  P  P  P  P  P  P  P 

boundary  conditions.  At  each  interface,  the  basic  laws  of  magnetostatics  (Equa¬ 
tions  (3a))  reduce  to  boundary  conditions  on  B  and  H  that  can  be  used  to  evaluate 
these  eight  constants.  The  normal  component  of  B  across  each  boundary  must  be 
continuous,  i.e.,  •  n^2  =  0  where  the  quantity  n^  is  the  unit  outward  nor¬ 

mal  to  the  surface.  This  provides  the  following  boundary  conditions  which  must  be 
satisfied  by  the  solutions  given  in  Equation  (160)  for  each  region 


B0I 

=  Bnii 

at 

zJ 

II 

Z3 

h* 

(161a) 

Bnii 

=  Bnm 

at 

n  =  n2 

(161b) 

Bnm 

Bniv 

at 

n  =  n3 

(161c) 

Bniv 

=  Bnv 

at 

Z3 

II 

ZJ 

(161d) 

However,  because  the  vector  potentials  in  each  region  are  functions  of 
P1(v),  we  can  simplify  Equation  (161)  to  constraints  on  A^  at  the  interfaces 


A.  T 

H 

=  Vi 

at 

n  =  n1 

(162a) 

=  %III 

at 

n  *  n2 

(162b) 

> 

M 

II 

at 

n  =  n3 

(162c) 

\iv 

=  \v 

at 

n  =  nA 

(162d) 
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The  second  set  of  boundary  conditions  states  that  the  theta  (8)  or  tangential, 
component  of  H  across  each  boundary  must  satisfy  the  relationship 


n12  X  (H2“H1)  =  Js 


(163) 


where  j^^which  equals  J^(8)^  is  the  real  surface  current  density  in  the  limit  of 
vanishing  width  between  the  two  regions.  Using  the  relationship  B  =  pH,  Equation 
(163)  can  be  expressed  as 


B02  B01 


=  3,1,  (9) 


(164) 


The  current  must  be  expanded  in  a  series  of  associated  Legendre  functions 

PL(v)  as  in  Reference  9.  The  form  of  the  current  is 
P 


OO 

J  Z  Vp  (cos  0) 

ve)  =  - 2 — ~ 

'  a(sinh  ri  +  sin  0)  2 


(165) 


where,  using  £  =  cosh  ri  and  v  =  cos  0,  V  can  be  shown  to  be 

P 


-(2p+l)  a 
2p(p+l) 


(^2-v2)^  pj(v)dv 


(166) 


For  the  two  current  bands  of  interest,  we  have 


J1  Z  V  pj(v) 

j  (0)  = - ELil - - - 

c2  2A 

a(  Sf-v  ) 


(167a) 
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where 


(167b) 


V  = 
P 


- ( 2  p+1 )  a 
2~p(p+l) 


P1(v)dv 

P 


and 


ve) 


J2  X  UnP^(v) 


ri  p  p 


a(^-v2) 


where 


U  = 
P 


~(2p+l)  a 
2p(p+l) 


and  =  cos  0^  and  =  cos  6^ . 


Referring  to  the  curl  in  Equation  (120),  we  can  write  Bq  in  the  form 


Be  55  (VxVe 


1  5(e3Aip)  _ 


e!e3  3n 


kW) 


h  9? 


\e-i) 


h  -i 


From  Equations  (164)  and  (169),  the  tangential  components  of  B  in  regions  ] 
V  must  satisfy  the  relationship 
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(168a) 


(168b) 


(169) 


through 


3E 

K2'1^  (  Vp1(?>  +  GpQJ®)  Pp1(v) 
\  / 

€- 

/  \ 

\  1  9 

35 

<52"1)!if”ppp1<S)  +  epqp1<5:i)pp1(v) 

r- 

3 


t 


yj<54> 


GpQp<V 


pJ(v)  = 


H  Q1(C, )pl(v) 
p^p  4  pv 


/  , 

/  V 

i  1  i  1  8 

\  Ul)attl-v2)h  8? 

(CU^KqJtoW) 

f= 

1 _  3 


yl/a(^2-v2)^  ^ 


(c'-i^fFpPjtt}  +  GpQJ(OjpJ(v) 


Jp2 (6) 


By  making  the  following  substitutions 


PP<E>  -^(V-dVcs)] 

<«>  ■tef'?2-1)4  <£«>] 
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( 17  If ) 


(17 lg) 


(171h) 


(172a) 


(172b) 


and  performing  algebraic  manipulation,  the  eight  boundary  conditions  can  be  simpli- 
f  ied  to 


Vp«l>  '  Vp'5!1  +  VpK1> 


(173a) 


B  PiJ(f 
P  p  l 


>  +  cdqd(5i)  1 + As,)] 

pp  \pi/Lpp  1J 


Jpl(e)  a (£^-v2)' 


pj(v) 


Vp<(2>  +  vJ({2)  ‘  Vj«2» +  y)hy 


(ti)  [v?«2>  +  ,.] 


\  i 


B  P^A,,)  +  C 

p  p  2 


Wf2>'_ 


+  Cn^<rn)  =  F  P2(L)  +  G  Q!(-  ' 
P  P  5  P  P  j  p  p  3  P  p  i 


,  ,  .  u  p'V,,)  +  e  QA(f_) 

h  >  J  P  P  3  p  P  3' 


]■  Ms 


.A,-  ,  .  A. 


F  P  '  (C  o )  +  G 
p  p  3  p  'p  3 


(173b) 


(173c) 


(173d) 


a  7  3c) 


( 17  3f ) 


f„p,!<g>  +  <=  q'«o1  - «.«!«.) 


(173,:) 


Because  the  eight  coefficients  can  be  determined  for  a  specified  problem  from 
Equation  (174),  the  potentials  Aj,  An>  Aj  ,  AJV,  and  Ay  in  regions  I  through  V 
can  be  completely  determined.  The  normal  (B^)  and  tangential  ( Bq )  components  of 
the  magnetic  induction  in  each  region,  I  through  V,  can  be  determined  by  using 
Equations  (120)  and  (124)  to  be: 


bot  =  -  L 


1  3 


J01  ^  77  ^  ^ 

P-l  a(^-v2) 


(175a) 


B  =  _  v  i  i_ 

0 1  ^  /  -  n\h  3v 


P=1  a (?2~V2) 


(i.vf  (Vp(5)pp<v> 


(175b) 


1  3 


Byi1  ^  /  2"  2\V* 
P=1  a(c2-v2) 


(42-l)  (BppJ(C)  +  cpqh5))  pV) 


(175c) 


Bnii  £7 


l  3 


(i-v2)  (B/J(5)  +  cpqj(0)  pJ(v> 


(175d) 


Beni  = k  H  (DpPj<?>  +  v’(0) Pp 


(v) 


(175e) 


BnIII  “  '£> 


p=l 


JL _  _3_ 

3v 


(1-p2)'  (yj‘s>  +  bpqJ<£))pJ<v) 


(175f) 
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(175g) 


(17  5h) 


(1751) 


(175j) 


APPENDIX  A 


FERROMAGNETIC  SPHERICAL  BODIES  IN  A  CONSTANT 
EXTERNAL  INDUCING  FIELD 

INTRODUCTION 

In  previous  work  Brown  and  Baker'7  derived  the  closed  form  mathematical  expres¬ 
sions  for  the  magnetic  flux  density  for  two  configurations  of  a  magnetic  spherical 
body  surrounded  by  a  stationary  current  band  of  azimuthal  spherical  symmetry.  The 
first  case  treated  was  for  an  infinitesimally  thin  stationary  current  band  surround¬ 
ing  a  spherical  magnetic  shell.  The  second  case  is  for  a  stationary  current  band  of 
finite  width  surrounding  a  solid  magnetic  sphere.  The  magnetic  bodies  were  assumed 
to  be  linear  and  homogeneous. 

The  problem  ot  the  magnetic  induction  for  an  infinitesimally  thin  current  band 
surrounding  a  spherical  shell  can  be  generalized  to  include  an  external  magnetic 
field.  The  superposition  principle  discussed  in  the  text  of  this  report  can  be 
used  in  these  two  cases  to  include  a  constant  external  magnetic  field.  The  magnetic 
induction  in  each  region  for  a  three-dimensional  magnetic  spherical  shell  in  an 
arbitrary  external  magnetic  field  is  added  to  the  magnetic  induction  for  the 
corresponding  region  for  the  spherical  shell  surrounded  by  and/or  surrounding  a 
stationary  current  band.  The  problem  of  deriving  the  magnetic  induction  for  a  cur¬ 
rent  band  of  finite  width  surrounding  a  solid  ferromagnetic  sphere  can  also  be 
generalized  to  include  an  external  magnetic  field  H  in  a  similar  manner.  Thus,  the 
magnetic  induction  for  a  ferromagnetic  spherical  body  in  an  external  magnetic  field 
■  ■.list  in  determined. 

both  constant  external  field  problems  were  solved  by  Nixon*1  of  the  ('.enter. 

’■'he  cloned  form  mathematical  solutions  for  the  magnetic  induction  for  both  constant 
external  field  problems  were  presented  in  Reference  6  in  Cartesian  coordinates. 
was  necessary  to  convert  these  mathematical  expressions  to  spherical  coordinates  to 
he  compatible  with  this  work. 

SOI. ID  FERROMAGNETIC  SPHERE  IN  AN  EXTERNAL  INDUCING  FIELD 

I  he  solid  ferromagnetic  sphere  in  a  constant  external  inducing  field  is  shown 

m  I  O’.urc  A.!.  The  permeability  of  the  solid  sphere  is  p  and  the  radius  of  the 

xpltiuu  iu  K  |  .  The  permeability  D  of  vacuum  that  is  external  to  the  sphere  is 

denoted  by  ..  .  The  constant  arbitrary  magnetic  field  is  designated  as  H 

1  "  o . 
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Figure  A. I  -  Ferromagnetic  Spherical  Solid  in  a  Constant 
External  Magnetic  Field 
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is  constant  in  the 


It  is  assumed  that  p ^  in  the  sphere  is  constant,  and  that  jj 
region  external  to  the  sphere.  Because  there  are  no  currents  in  any  region  in  the 
problem,  the  magnetic  field  H  can  be  expressed  as  the  negative  of  the  gradient  of  a 
magnetic  scalar  potential  tf/  in  regions  I  and  II,  respectively. 

Ht  =  -  V <t>  for  0  <  r  <  R.  (A. la) 

1  Im  —  —  i 


where 


H  =  -  71  for  R-.  <  r  < 

II  I Im  1  — 


(A. lb) 


BI  =;;2HI 


(A.lc) 


BII  *  "lHII 


(A. id) 


The  major  step  toward  solving  this  problem  is  to  determine  the  solutions  of  the 
scalar  Laplace's  equation  in  regions  I  and  II  which  satisfy  the  boundary  conditions 
at  r  =  R^.  In  terms  of  B  and  H  the  magnetostatic  boundary  conditions  are 


(B  -B  )  •  n19  =0  at  r  =  R 


(A. 2a) 


n^  x  C  H i  — H i )  =  0  at  r  =  Rj 


(A. 2b) 


where  n^?  is  the  unit  vector  normal  to  the  surface  of  the  sphere  from  region  I  to 
II . 

The  general  expression  of  <J)  and  <b  which  satisfy  Laplace's  equation  in 

Im  Ilm 

regions  I  and  II  are: 


*Im 


Ar  sin  h  cos  +  Br  sin  0  sin  +  C  r  cos  0 


for  0  <  r  ^  Kl 


(A. 3a) 
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Ilm 


-2  -2  -2 
=  Dr  sin  8  cos  +  Er  sin  8  sin  +  Fr  cos  6 


-  H  r  sin  8cos  tii  -  H  r  sin  8  sin  tf;  -  H  r  cos  8 
ox  T  oy  r  oz 


for  _<  r  <  “ 


(A. 3b) 


The  coefficients  are  determined  by  the  magnetostatic  boundary  conditions  (see 
Equations  (A. 2a)  and  (A. 2b)  and  are: 


A  =  - 


VA 


(A. 4a) 


B  =  - 


3VJlHoy 

U2+2Ui 


(A. 4b) 


C  =  - 


3M  -j  H  , 
1  oz 

U2+ 2b1 


(A. 4c) 


D  = 


R  (y  -m.)H 
1  2  1  ox 

(U2+2y1) 


(A.4d) 


E  = 


R1^2~M1^H- 


22 


(^+2^) 


(A.4e) 
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(A.4f ) 


F  = 


R1(1J2“M1)H 

~^T 


oz 


For  details  of  this  derivation,  the  reader  should  consult  the  work  of  Nixon. ^ 
The  authors  have  transformed  his  mathematical  expressions  for  the  magnetic  induc¬ 
tion  in  regions  I  and  II  from  Cartesian  to  spherical  coordinates.  This  makes  the 
expressions  for  the  magnetic  induction  compatible  with  the  present  work  of  this 
report.  These  mathematical  expressions  are  in  regions  1  and  II; 


B^(r,9,4')  =  -  M  2^r  ^  sin  ®  C0S  'f*  "*■  B  sin  9  sin  *  c  cos 


-  U2eg(A  cos  0  cos  tf/  +  B  cos  0  sin  ip  -  C  sin  0) 


-  J-i .(~A  sin  ^  +  B  cos  i|0  for  0  r  £ 


(A. 5a) 


h^e^(2Dr  ^sin  0  cos  iji  +  2  Er  ^ 


sin 


_3 

sin  4*  +  2  Fr  cos  6 


+  H  sin  o  cos  \b  +  H  sin  0  sin  ui  +  H  cos  0) 
ox  r  oy  oz  ' 


/\  _3  _3  _3 

-  p  e^(Dr  cos  0  cos  +  Er  cos  0  sin  <Jj  -  Fr  sin  0 
Note:  Above  equation  continued  on  next  page. 
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-  H  cos  0  cos  ti  -  H  cos  0  sin  <1  +  H  sin  0) 
ox  oy  oz 


^iV-Dr' 


sin  +  Er  cos  ^  4* 


H  sin  ip  -  H  cos  <10 
ox  oy 


(A. 5b) 


where 


B  =  -  pV<t> 

m 


(A. 5c) 


and 


34> 

_  /x  m 
"  er  9r 


/\ 

e0 


3  4> 

m 
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d<P 


r  s  in  0  dip 


(A.5d) 


FERROMAGNETIC  SPHERICAL  SHELL  IN  AN  EXTERNAL  INDUCING  FIELD 

The  problem  of  the  spherical  shell  is  similar  to  the  problem  of  the  solid 
sphere.  The  inner  radius  of  the  spherical  shell  is  R^  and  the  outer  radius  is  R£ 
(see  Figure  A. 2).  The  permeability  of  the  magnetic  material  in  the  shell  is  M2 
and  the  permeability  UQ  of  vacuum  that  is  internal  and  external  to  the  shell  is 
denoted  by  M-^ •  The  constant  external  magnetic  field  is  designated  by  Hq. 

The  problem  of  deriving  the  closed  form  mathematical  expressions  for  the 
magnetic  flux  density  in  each  of  the  t’lree  regions  (I  through  III)  was  worked  out 
in  detail  by  Nixon.  The  problem  was  solved  in  a  method  exactly  analogous  to  the 
problem  of  a  solid  sphere  in  a  constant  external  magnetic  field.  For  details  of 
the  derivation  consult  Reference  6. 

The  general  expressions  for  the  magnetic  scalar  potential  4>m  in  regions 
I  through  III  are; 

<b  =  Ar  sin  0  cos  ip  +  Br  sin  0  sin  ip  +  Cr  cos  0  for  0  <  r  <  R.  (A.  6a) 

Im  —  —  i 


*TIm 


-2  -2 
=  (Gr+Hr  )  sin  0  cos  tj>  +  (Ir+Jr  )  sin  0  sin  <j) 


_? 

+  (Kr+Lr  )  cos  0  for  _<  r  _<  R^ 


(A. 6b) 
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L23 


Illm 


-2  -2  -2 
=  Dr  sin  8  cos  i|)  +  Er  sin  0  sin  ^  +  Fr  cos  0 


-  H  r  sin  0  cos  vli  —  H  r  sin  0  sin  ^ 
ox  oy  r 


-  r  cos  0  for  <  r  <  00  (A.  6c) 


The  coefficients  are  determined  by  the  usual  magnetostatic  boundary  conditions 
on  the  spherical  surfaces  at  r  =  R^  and  r  =  R2.  The  coefficients  determined  by 
this  method  are: 


H 


!^1 


1-1 


(2y i+M  2 ) (2p2+p1) 


-u2) 


(U2-M1)R1 


ox 


(A. 7a) 


(2y  +y  )  H 

G  =  - -  ---  -  (A. 7b) 

(y  -y  )rJ 


D  =  GR^  +  H  +  H  R? 
2  ox  2 


(A. 7c) 


A  =  G  +  HR, 


(A.7d) 


J  =  -3y. 


(2y1+y2) (2y2+yL) 


2(yL-y2) 


(y2-M1)R1 


oy 


(A. 7  e) 


124 


1 


J 


(A . 7f ) 


(2p2-HJ1) 

(y2-u1)  R 


3 

1 


E  = 


+  J  +  H 

oy 


R 


3 

2 


B 


I  +  JR 


-3 

1 


L  =  -3y. 


(2y j+u9) (2u2+m1) 


(y0-u,)  R, 


(2m1-M2) 


-1 


K  = 


(2y2+P1) 

(VJ2~tJl)  R] 


3  3 

—  KR,,  +  L  +  H  R.. 
2  oz  2 


C  =  K  +  I.R  3 


The  authors  have  transformed  Nixon's*3  mathematical  expressions  for  the 
induction  in  regions  1  through  Ill  from  Cartesian  to  spherical  coordinates, 
m.-ike:;  tin-  expressions  for  the  magnetic  induction  compatible  with  the  present 
this  report.  These  mathematical  expressions  in  regions  1  through  111  are: 


125 


-  - - *  —  . mJti, 


(A.7g) 


( A .  7  h ) 


( A .  7  i ) 


(A. 7 j ) 


(  A .  7  k ) 


(A. 7;  ) 

magnet ic 
This 
work  of 


Bj(r,0,^)  =  -  M^e^(A  sin  0  cos  ip  +  B  sin  0  cos  >p  +  C  cos  0) 


-  y/eX(A  cos  0  cos  to  +  B  cos  0  cos  -  C  sin  0) 
i  0 


-  y 


(-  A  sin  ijj  +  B  cos  ip) 


for  0  r  £ 


B^r.e,*)  = 


/\ 
-  u2e 


_3 

2r  H)  sin  0  cos  ip 


+  (I-2r  3J)  sin  0  sin  ^  +  (K-2r  3L)  cos  0 


-  y2'e^)  J^(G+Hr  3)  cos  0  cos  ip  +  (I+Jr  3) 


(cos  0  sin  (|0 


(K+Lr-3) 


sin 


U  ^ 
M2 ip 


-(G+Hr  3)  sin  ip  +  (I+Jr  3) 


for  R  _<  r  £  R2 


(A. 8a) 


(A. 8b) 
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|2Dr  ^  sin  0  cos  1])  +  2Er  ^  sin  0  sin  iJj 


+  2Fr  cos  8  +  H  sin  0  cos  ip  +  H  sin  0  sin  \p  +  H  cos  0 
ox  oy  oz 


/\  y.  “3  «  Q 

^le0  r  cos  0  cos  \p  +  Er  cos  0  sin  ip  -Ft  sin  0 


H  cos  0  cos  ip  -  H  cos  0  sin  ij)  +  H  sin  6 
ox  oy  oz  J 


-  y,’e',  f-  Dr  ^  sin  il)  +  Er~^  cos  ip  +  H  sin  ip  -  H  cos  4M  (A. 8c) 
1  ^  1_  ox  T  oy  J 


where 


B  =  -  yV4> 


(A.8d) 


3<t>  'el  9$  'e',  34> 

v  — E  +  JL  _JL  +  _ !fe _  _J 

r  3r  r  30  r  sin  0  3^ 


(A.8e) 


APPENDIX  B 


DERIVATION  OF  THE  COEFFICIENTS  OF  THE  VECTOR  POTENTIAL 
FOR  A  SOLID  SPHERE  SURROUNDED  BY  AN  INFINITESIMALLY 
THIN  SPHERICAL  CURRENT  BAND  AND  REDUCTION  OF  THE 
MAGNETIC  VECTOR  POTENTIAL  TO  THAT  OF  A  THIN 
BAND  IN  VACUUM  WHEN  IN  THE  LIMIT  P2  EQUALS  l-^ 

DERIVATION  OF  THE  COEFFICIENTS 

In  this  appendix  the  coefficients  are  derived  for  the  vector  potential  in 
regions  I  through  III  for  a  ferromagnetic  sphere  surrounded  by  an  infinitely  thin 
current  band.  For  a  detailed  discussion  of  this  ferromagnetic  problem,  see  the 
section  in  t  lie  text  of  the  report  entitled  "Solid  Sphere  Surrounded  by  an 
Infinitesimally  Thin  Spherical  Current  Band".  The  magnetic  vector  potential  in 
each  region  is  given  by; 


AT  =  A 


i  I 


z  (■ <cos  B> 

P=i 


(B.la) 


A  =  A. TT 
II  <MI 


L 

p=i 


A  „rp  + 


JlL 


p2  (p+1) 


P  (cos  6) 
P 


(B. lb) 


AI 1  I  A',’TII  2Z 

P=1 


J* 


.(p+1) 


P  (cos  0) 
P 


(B.lc) 


The  coefficients  (A  . 

Pi 

substituting  these  eci 


)  and  (B  in  Equations  (B.la)  through  (B.lc)  are  obtained  by 


nations  into  tlie  boundary  conditions  (Equations  (B.2a)  through 


(B.2c) 


-  v;  r  If  (rAu}  +  ^  7  h  (rV  =  0  at ! =  Ri 

7  97  (rAIII)  +  7  a7  *rAII*  =  J(0)  at  r  =  R2  (B.2d) 


After  appropriate  substitutions  of  Equations  (B.la)  through  (B.lc)  into  Equa¬ 
tions  (B.2a)  through  (B.2d),  the  following  boundary  value  equations  are  obtained. 


'pirip  -  [  Vip  +  BP2Er(p+1>] 


(B.3a) 


[Ap2R2P  +  Bp2R2_<P+1)]  ‘  Bp3R; 


-(P+1) 


(B .3b) 


0  (B.3c) 


[k  -  PBp2R2-lP+2)J-  £ 


J  (0) 


(B.3d) 


P  (cos  0) 
P 


These  algebraic  equations  provide  four  simultaneous  equations,  with  four 
unknowns,  which  can  be  solved  for  the  coefficients  A  ^  and  B  ^  by  algebraic  mani¬ 
pulation  . 

Solving  Equation  (B.3a)  for  A  and  Equation  (B.3b)  for  B  .  we  have 

pi  po 
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Apl  Ap2  +  Bp2Rl 


-(2p+l) 


(B.4) 


A  RP  +  B  ,R:(p+1) 

B  .  =  -Eli _ -J2  L 


p3 


R, 


-(p+1) 


(B .  5) 


Solving  Equation  (B.3c)  for  A 

P1 


gives 


R-  (2  p+1)  -| 
Bp2pRl 

(p+1) 


(B .  6) 


which, 


when  substituted  into  Equation  (B.3a)  and  solved 


for 


yields 


(B.7) 


Substituting  Equations  (B.5)  and  (B.7)  into  EquaL  on  (B.3d)  gives  the  expression 
for  B 

P2 
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Simplifying  Equation  (B.7)  by  using  Equation  (B.8)  yields 


(0) 

-L  D 


p2  (2p+l)  R^P_1)p1(cos  6) 

2  p 


(B.  9) 


The  following  expression  is  obtained  for  A  ^  from  Equation  (B.4)  after  substituting 
Equations  (B.8)  and  (B.9)  for  B  and  A^2,  respectively,  as 


U.J  (0) 

^  =  _ LP - 

pl  (2P+l)R^P_1)Pp(cos  6) 


l-*I 


y2\VlJp«> 

P1  /  Bp(cos  6) 


(P-1)/  / 

(2p+l)R2  I  1  +1 


(P+D 


(B .10) 


The  mathematical  solution  for  B  ^  in  terms  of  known  quantities  is  obtained  from 
Equation  (B.5)  by  substituting  the  previously  obtained  expressions  for  A^2  (Equa¬ 
tion  (B.9))  and  for  B  (Equation  (B.8)). 

p2 


i  J  (0)R$2p+1) 
-LP_.  .2 _ 

P^(cos  0) 


1  _  M2\/Vp(9)  ^ 

Wl/\p\cos  6), 


(2p+l)R 


(p-1) 


(  R-(2p+1)  )  (2P+l)R^p-1)  (l+^ 


(P+1) 


(B. 11) 
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After  the  numerical  value  for  B  ^  is  calculated  on  the  computer  for  a  specific 

P3 

problem,  the  numerical  values  for  the  other  coefficients  can  be  obtained  from  the 
following  equations: 


M  mijp(8) 


U1 )  P 1 ( co s  6) 


P3 


(-R  (2p+1)^  (2p+l)R^p_1) 


(*•© 


p+1 


(B . 12a) 


(6) 

_ L_p_ _ 

P^cos  0) 

_ 2_ 


P  2 


(2p+l)R, 


(P-D 


(B12 . b) 


L  _  piJP(Q) 


Pi 


11 1 )  P1  (cos  B) 


(Pp+IJR^  ^^Pp(cos  9) 


(2p+l)RP 


rLfi 

-  y  /  p+i  / 


(B  ,12c) 


REDUCTION  OF  THE  POTENTIALS  WHEN  p  EQUALS  p] 

The  vector  potentials  for  this  problem  should  reduce  to  those  of  an 
infinitesimally  thin  current  band  in  a  homogeneous  medium  with  permeability  p^ 
in  the  limit  as  p  =  p  (see  Figure  B.l).  In  this  limit  the  coefficients  should 
assume  the  following  form 
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A=A,;B=0;B,*0 
pi  p2‘  p2  p3 


(B . 13) 


where  A^  and  B  ^  should  reduce  to  the  coefficients  for  the  potentials  in  the  two 
regions  for  the  spherical  band  problem  of  Appendix  B  in  Reference  7.  One 
immediately  observes,  from  Equation  (B.12a)  that  B  ^  =  0  when  =  From  Equa¬ 

tions  (B . 12b)  and  (B.12c)  we  see  that 


r  u,j  (9) 

_ L_e _ 

P^(cos  0) 


(2p+1)rJp 


(B . 14a) 


A  .  =  A  _ 
pi  P2 


(B . 14b) 


From  Equation  (B.ll),  in  the  limit  of  =  h-j.'  we  have 


P,J  (9) 

_ L_e _ 

R2  ^P+^\2p+l) (cos  0) 


(B . 15) 


Rewriting  A  .  we  have 
Pi 


A  =  B  R  “(2p+1) 
pi  P3  2 


(B .16) 
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This  means  that  in  the  three  regions  the  components  of  the  vector  potentials 
used  in  Equations  (B.l): 


£  (APirP)  pp(cos  6) 

p=i 


Vi  =  £  VP+  ifriy  pJ(cos  e) 

p=l  L  r  J 


Vn  "  ^  r(p+D  pP(cos  e) 

P=1  L  r  J 


reduce,  when  =  to  t*ie  f°rm 


oo 

&(■ 


'y2^i 


rP  P^Ccos  0) 
P 


=  £  /B 


— 1 

|Bp3  jr_(p+1)  pJ(cos  6) 

\  y2=-yl '  J 


(B . 17a) 


(B.17b) 


(B .  17c) 


(B .18a) 


(B . 18b) 


The  mathematical  expressions  for  A  ^  and  B  ^  (see  Equations  (.B.16)  and  B.15), 
respectively),  for  the  ferromagnetic  sphere  surrounded  by  a  thin  current  band  in 
the  limit  as  =  P-j.*  are  t^e  same  as  f°r  t^ie  coefficients  A  ^  and  B^  (see  Refer¬ 
ence  7,  Appendix  B) ,  respectively,  for  the  components  of  the  vector  potentials  in 
the  regions  of  the  current  band  in  vacuum.  For  comparison,  the  coefficients  for 
the  current  band  problem  are: 


13  6 


A 


pl 


B  „r7(2P+1> 


u,j  (0) 

B  =  - i-E- - 

p2  pJ(cos  6)R^(p+2)(2p+l) 


and  the  coefficients  for  the  ferromagnetic  sphere  problem  with 


=  B  ®-<2P+D 
p3  2 


y,J  O) 

B  =  - - 

p3  P3(cos  0)R~(p+2)(2p+l) 


It  is  noted  when  making  the  comparison,  R2  must  be  set  equal  to  R^ . 


(,B.19a) 


(B . 19b) 


M1  are 


(B.20a) 


(B. 20b) 
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APPENDIX  C 


DERIVATION  OF  THE  COEFFICIENTS  OF  THE  VECTOR  POTENTIAL 
FOR  AN  INFINITESIMALLY  THIN  CURRENT  BAND  SURROUNDED 
BY  A  FERROMAGNETIC  SPHERICAL  SHELL  AND  THE  REDUCTION 
OF  THE  MAGNETIC  VECTOR  POTENTIAL  TO  THAT  OF  A  THIN 
BAND  IN  VACUUM  WHEN  IN  THE  LIMIT  p2  EQUALS  ^ 

DERIVATION  OF  THE  COEFFICIENTS 

In  this  appendix  the  coefficients  are  derived  for  the  vector  potential  in 
regions  I  through  IV  for  a  ferromagnetic  sphere  surrounding  an  infinitely  thin  cur¬ 
rent  band.  For  a  detailed  discussion  of  this  ferromagnetic  problem,  see  the  sec¬ 
tion  in  the  text  of  the  report  entitled,  "Hollow  Sphere  Surrounding  an  Infinitesi¬ 
mally  Thin  Spherical  Current  Band".  The  magnetic  vector  potential  in  each  region 
is  given  by: 


Ax  =  A,  T 
I  ^1 


=  L 


p=i 


A  r* 
Pi 


P  (  cos  6) 
P 


(C. la) 


AII  %I1 


z 

P=1 


|A  rP  +  — 1 £ 


p2‘  '  r(p+l) 


P  (cos  6) 
P 


(C.lb) 


AIII  =  Avm 


L 

P=i 


A  rP+^ 


p3  r(p+D 


P  (cos  0) 
p 


(C.lc) 


A  =  A 
IV  iJjIV 


L 

P=i 


j± 


.(p+i) 


P^(cos  0) 
P 


(C.ld) 


Ihe  coefficients  (A  ,  and  B  ,)  in 

pi  pi 

substituting  these  equations  into 
t  hrough  ( C . 2 1 ) ) . 


Equations  (C.la  tlirough  C.ld)  are  obtained 
the  boundary  conditions  (Equations  (C.2a) 


by 
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PRECEDING  PAGE  BLANK-NOT  FILMED 


(C.2a) 


AI  AII  r  "  R1 


AII  AIII  r  =  R2 


AIII  AIV  r  R3 


k  7  h  <rAn> + 1  ?  If  <rV  •  J<0> 


r  =  R, 


V2  r  3r  (rAIll)  +  y  r  3r  (rAIiJ  ~  0 


r  =  R„ 


—  —  —  (rAT,,)  +  — ■  —  -r—  (rATTT)  =  0  r  =  R„ 

r  3r  IV  y2  r  3r  v  III'  3 


After  appropriate  substitutions  of  Equations  (C.la)  through  (C.ld) 
tions  (C.2a)  through  (C.2 f),  the  following  boundary  value  equations  are 


*plRl  '  Ap2Ef  +  Bp2RIlP+1) 


\A  +  Bp2R;l,,+1>  •  Ap3R2  +  Bp3R2<P+1> 


A  RP  +  B  r:(p+1)  =  B  ,R“(p+1) 
p3  3  p3  3  p4  3 


(C.2b) 

(C.2c) 


(C.2d) 


(C.2e) 


(C.2f) 


into  Equa- 
obtained . 

(C.3a) 

(C. 3b) 

(C.3c) 
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A  (p+l)R 


(P-1) 

1 


1 


(P+1)rJP_1) 


P^(cos  j ) 


(C .  3d) 


Ap3(p«)R<f-1) 


l’Bp3R)(P+2) 


Ap^ (p+1) R 


(p-1) 

2 


PE 


p2 


-  (p+2) 


=  0 

(C . 3e) 


■>V)(p+2) 


AP3(P+1)R3P  1} 


»BP3R3(P+2) 


=  0 


(C.3f) 


The  algebraic  equations  provide  six  simultaneous  equations  with  six  unknowns  to 
which  can  be  solved  for  the  coefficients  A  and  B^,  by  algebraic  manipulation. 
Solving  Equation  (C.3a)  and  Equation  (C.3d)  for  A  respectively,  we  have 


B  ?R-(2P+1> 
P2  1 


(C.4a) 


p  J  (0) 
-1  -£ _ _ 


P  (cos  0) (p+1) 


R"(p_1)+  A  2 
1  P2 


*  p^(Zp+1> 

S- I - — - 

(p+1) 


(C ,4b) 


Equating  (C.4a)  and 


(C.4b)  and  solving  for  B 


P2 


yields 


M  J  (0)RP  2 

B  =  - -  =  j'(B)  (C.5) 

(2p+l)P  (cos  0)  P 


141 


f 


Using  Equation  (C.3c)  to  solve  for  B  .  we  have 

p4 


Jp4‘ 


A  ,K«P+1> 
p3  3 


+  B 


p3 


(C.6) 


and  similarly,  using  Equation  (C.3f)  to  solve  for  B 


p4 


we  have 


(C.7) 


Now,  upon  equating  Equations  (C.6)  and  (C.7)  and  solving  for  B 


p3 


we  obtain 


(C.8) 


where 


Solving  Equation  (C.3b)  for  A 


p2 


and  substituting  Equation  (C.3)  for  B 

P2 


yields 


j'(0)r:(2p+1) 

P  2 


+  A 


p3 


+  B  _R"(2p+1) 

p3  2 


(C.9) 


Using  Equations  (C.5),  (C.8),  and  (C.9)  in  Equation  (C.3e)  yields  the  following 

expression  for  A  . 

p3 
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1 


(p+l)R$P  X) 
l  i 


[x J(p+1) 


J^g)(2p+l)R"(p+2)j 


~(p+2)  _JL  (p+1)rCp-i> 
U2  1 


where 


jp(6) 


vP<««r 

(2p+l) (cos  0) 
P 


I  he  constants  have  now  been  found.  After  the  numerical  value  of  A  .  is 
r  .  p3 

for  a  specific  problem,  the  numerical  values  for  the  other  coefficients 

obtained  from  the  following  equations: 


B 


P2 


(0) 


v  -  vM 


A  2  =  -  d' (H)R?(2p+1)  +  A  +  B  R_(2p+1) 
p2  p  2  P3  p3 


*  ,  =  A 
pi  p2 


B  ?R-(2p+1) 

p2  1 


[xK<p+2>] 

(C.10) 


calculated 
can  be 


(C.lla) 


(C.llb) 


(C.llc) 


(C.lld) 


143 


V  '  ap3R32P+1>  +  Bp3  (C-lle) 


REDUCTION  OF  THE  POTENTIAL  WHEN  P2  EQUALS  y1 

When  y^  is  set  equal  to  y2  the  above  ferromagnetic  problem  reduces  to  that  of 
finding  the  potentials  in  the  two  regions  of  a  simple  current  band  (see  Figure 
(B.l),  because  the  ferromagnetic  shell  will  now  have  a  permeability  y^  equal  to  that 
of  the  homogeneous  medium  with  a  permeability  y^.  In  this  limit  the  coefficients 
should  assume  the  following  form: 


Apl  *  0  Bp2  =  V 


A  0  =  0  B  ,  =  B  , 

p2  p3  p4 


A  ,  =  0  B  ,  +  0 

p3  P^ 


(C.12) 


where  A^  and  B^  should  reduce  to  the  coefficients  for  the  potentials  in  the  two 
regions  for  the  spherical  band  problem  (see  Reference  7,  Appendix  B) .  If  the 
coefficients  assume  this  mathematical  form  it  will  prove  that  the  mathematical 
forms  of  the  coefficients  for  the  spherical  shell  surrounding  a  thin  current  band 
are  mathematically  correct. 

The  coefficient  A^  will  now  be  evaluated  when  the  limit  is  taken  with  y^  =  y^ 
which  causes  to  approach  infinity  (see  Equations  (C.10)). 


A  »  =  limit 

P3  r  n 


M2=M1 


±- /(2P+m;(r+2> 


^  <p+i)»2p-1h  L  [x](p+l)R2-^«-  i_(p+i,R(P-l)+  I_  p[x]R 


-(p+2) 

2 


(C.13) 
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0 


p2^l 


where 


J 

P 


M,J  (0) 

_ i_E _ 

(2p+l)R~(p+2V(cos  0) 


The  expression  for  B  „ 

p2 


(see  Equation  (C.lla))  is 


(C.14) 


The  expression  for  B 

M 


when 


p,  is  still  J  since  J  is  not  a  function  of 

1  p  p 


P2=Pl 


(C . 15) 


The  expression  for  B 

P3 


(see  Equation 


(C.llb)) 


is 


v  ■  vM 


(C.16) 
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where  y,,  =  can  be  written  as 


P3 


Mir  v2p+1>r"2(p+2) 


=  Jlimit 


^  H 


i-(p+l)R<P-  !>+  i—  [x|  (p+1)r;<P+2L  i-(p+1)E<l'-l)+  I 

M1  z 


&imit 


^-[x]/  (2p+l)R"(p+2) 


[x]->  -  (  p1[X]^2P+1')R2 


- (p+2) 


Jp3 


=  J 


U2=Ml 


The  expression  for 


(see  Equation  (C.llc) 


is 


i 

J  R 
P 


-(2p+l) 

2 


+  A 


P3 


+  B 


P3 


-(2p+l) 

2 


when  u-  =  y, ,A  „  can  be  written  as 


The  expression  for  A 


pi 


(see  Equation  (C.lld)) 


is 


+  B  ,R 
P2 


~(2p+l) 

1 


Pi 


=  A 


U2-lJl 


-(2p+l) 


\r-(2p+‘) 


u2=Ui  Vyi 


pA 


-(2p+l) 

1 


P2=Mi 


^1 


Pi 
y2=Ul 


(j'p) 


,-(2p+l) 


1A8 


This  means  that  in  the  four  regions,  the  components  of  the  vector  potentials 
used  in  Equation  (C.l) 


V  =  E  [VP]  PJ(C0S  0) 

P=1 


A  =  E  A  _rP  +  — pl(cos  9) 
P2  r(P+1)  P 

d=1  r 


A,  =  V]  A  rP  +  ■7Pxi"\  p  (cos  9) 
i^III  ^  p3  r(p+l)  p 


00  r  b  /  i  , 

A  .  y  p1-  p1 

Vv  pti  [r(p+DJ  P 


(cos  0) 


reduce,  when  ^  =  v^,  to  the  form 


A  =  E  |A  , 
<J»I  ,  Pp 
P=1  \ 


rP  P^(cos  0) 


U2=P1 


''I'll,  III,  IV 


E  B 


U2=Pi 


r  ^P+^  (cos  0) 
P 


(C.25a) 


(C ,25b) 


(C.25c) 


(C.26d) 


(C.26a) 


(C.26b) 


The  mathematical  expressions  for  A^  and  B  ^  (see  Equations  (C.24)  and  (C.23), 

respectively)  for  the  ferromagnetic  spherical  shell  surrounding  a  thin  current  band, 

in  the  limit  as  =  y, ,  are  the  same  as  the  coefficients  A  n  and  B  _  (see  Reference 
Z  1  pi  pZ 

7,  Appendix  B) ,  respectively,  for  the  components  of  the  vector  potentials  in  the 
regions  of  the  current  band  in  vacuum.  For  comparison,  the  coefficients  for  the 
current  band  problem  are 


B 


-(2p+l) 

p2  1 


(C.27a) 


y.J  (0) 

B  =  - i_£L - 

p2  Pp(cos  0)R"(p+2) (2p+l) 


(C.27b) 


and  the  coefficients  for  the  ferromagnetic  shell  problem  with  are 


,-(2p+l) 

1 


(C .28a) 


U,J  O) 

_ 1  P _ 

P1(cos  e)R7(p+2)(2p+l) 
pl 


(C.28b) 
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APPENDIX  D 


DERIVATION  OF  THE  COEFFICIENTS  OF  THE  VECTOR  POTENTIAL  FOR  A 
SPHERICAL  SHELL  WITH  INTERNAL  AND  EXTERNAL  INFINITESIMALLY 
THIN  SPHERICAL  CURRENT  BANDS  AND  THE  REDUCTION  OF  THE 
MAGNETIC  VECTOR  POTENTIAL  TO  THAT  OF  A  THIN  BAND 
IN  VACUUM  WHEN  IN  THE  LIMIT  p2  EQUALS 

REDUCTION  OF  THE  COEFFICIENTS 

In  this  appendix,  the  coefficients  are  derived  for  the  vector  potential  in 
regions  I  through  V  for  a  ferromagnetic  spherical  shell  with  internal  and  external, 
infinitely  thin,  spherical  current  bands.  For  a  detailed  discussion  of  the  ferro- 
magnetic  problem,  see  the  body  of  this  report.  The  components  of  the  magnetic  vec¬ 
tor  potential  in  each  region  are  given  by 


AT  =  A . t 

I  tpi 


£  fAPirP]pp(cos  e> 

p=l 


(D. la) 


AII  ~  Ai(jII 


£ 

P=i 


o 

P  +  _E2, 


A  _rr  +  . . 

p2  r(.P+l) 


P  (cos  0) 
p 


(D . lb) 


A  =  A 
III  >111 


£ 

p=i 


A  ,rP  + 


p3  (p+1) 


P  (cos  0) 
P 


(D.lc) 


Aiv  A4)iv 


=  A 


1|)V 


£ 

P=1L 


A  ,  rP  + 


je!l 


p4  (p+1) 


P  (cos  0) 
P 


£ 

p=i 


.(p+1) 


P  (cos  9) 
P 
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(D. Id) 


(D.le) 


A 


The  coefficients  in  Equations  (D.la)  through  (D.le)  are  obtained  by  substituting 
these  equations  into  boundary  condition  (Equations  (D.2a)  through  (D.2h). 


AI  =  AII  at  r  =  Rx 


AII  =  AIII  at  r  =  R2 


(D. 2a) 


(D. 2b) 


AIII  =  AIV  at  r  =  R3 


(D. 2c) 


Aiv  =  \  at  r  =  R4 


(D. 2d) 


^  7  17  <rAii>  +  57  7  h  <rAx>  '  Ve>  r  '  Ex  (D'2e) 


1  1  8 


1  1  a 


V  r  gr  ^rAIII^  +  r  ir 


(rAXI)  =  0 


r  =  R 


2  (D.2f) 


1  13/AN4.1  18,.  A  A 

—  —  tt-  (rA_„)  +  —  —  -s—  (rATTT)  =  0 

r  8r  IV  r  9r  111 


r  =  R 


3  (D.2g) 


^  7  h  “V +  77  7  If  <tAxv>  ■  J2(8>  r '  e4  (D-2h> 
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After  the  appropriate  substitutions  are  made,  the  following  equations  are  obtained: 


VR?  ■  [v< +  Vh(p+l)] 


P  A  n  n-(p+1)l  _  ("a  dP  r  D~(p+l)l 


A  ,R*  +  B  -R  =  [a  „R‘;  +  B  „R 

p2  2  P2  2  J  L 


p3  2  p3  2  J 


Ap3R3  +  Bp3E3<P+1)]  *  [Ap4R3  +  Bp4R3<P+1)] 


Ap4R4  +  Bp4R4"<P+1>]  ’  [VR3 


-  [VR4-<P+1>] 


Hv<p+1,s!p~1)  '  p0p2R3"<p+2)]+ ^[<p+1> ApiR{ 


(p-1) 


V6> 

P^Ccos  6) 
P 


1  rA  fp+DR^P-^  -  pB  JC(P+2)1+  ’  f a  _  „R  o-(P+2) 


1  > 


pi'1 - 2  “  '’VV  ^  [A,,2 


(ptl)R^  -  -  Pi^,3, 


Ap4(p+l>R5'"l>  -  pBp4Rj<P+2>l+  ~ 


A„3(p-U)I!<p-1)  -  pnp3R;(p+2) 


■  (p+2 )r 


PR 

p  J 


V  - 

J  Up 


V(|,tl)R4M)  '  'V^2> 


J  ,,(u) 

P2 _ 

P  V  os  P) 


(D.3a) 

(D .3b) 


(D.3c) 


(D . 3d) 

(D.3e) 


=  0 
(D.3f) 

0 

(D.3g) 

(D.3h) 
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These  algebraic  equations  provide  eight  simultaneous  equations  with  eight  unknowns, 
which  can  be  solved  for  the  coefficients  and  B  ^  by  algebraic  manipulation. 

Solving  Equations  (D.3h)  and  (D.3d)  for  B  ^  and  equating  the  results  to  solve 
for  ap4»  yields 


UlJp2(9) 


p4  (2p+1)rJp“1)p1(cos  0) 

4  D 


5  H 


(D.4) 


Solving  Equations  (D.3b)  and  (D.3f)  for  B 
Ap2 »  yields 


p3 


and  equating  the  results  to  solve  for 


(D.5) 


Solving  Equation  (D.3a)  and 
B^2>  yields 


(D.3e) 


for  A 


Pi 


and  equating  the  results  to  solve  for 


p2 


(2p+l)R 


~(p+2) _1 


P~(cos  0) 
P 


■H 


ID.  6) 


Solving  Equation  (D.3c)  and  (D.3g)  for  B  ^ 

AP3’  yields 


and  equating  the  results  to  solve  for 


V>  ■  “p3  ["]  +  MM  =  M 
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p  1+  3  K 

V  P  / 

L.  V  ' 


H  ■  v 

Now,  using  Equation  (D.3b)  and  the  results  of  Equations  (D.4)  through  (D.9)  to 
solve  for  B  yields 

,  tfHH  -  Rz M [s] M  -  MM 
P3  E2  M[t]  -  R2<P+1>  -  R2  [“] 


(D. 14) 


(D.15) 


Now  B 


P4» 


and  A  ,  can  be  found  from 
pl 


p4  p3 


,(2p+l) 

3 


+  B 


p3  ”  p4 


,(2p+l) 

3 


(D. 16) 


(2p+l) 

'4 


+  B 


p4 


Vp2(6) 


P1(.cos  0) pR.  (p+2) 
P  4 


(D.17) 


^  ■,  =  A  - 
pl  p2 


b  „r7C2p+1) 

p2  1 


(D.18) 


The  constants  have  now  been  found.  After  the  numerical  values  of  B  A  ,  and  B 

p3  p4  p2 

are  calculated  for  a  specific  problem,  the  numerical  values  for  the  other  coeffi¬ 
cients  can  be  obtained  from  the  following  equations 
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(  D .  1 9a ) 


flP3  -  v  M +  M  H  2  H 


P3 


(2p+l) 

3 


+  B 


P3 


(2p+l) 

3 


(D.19b) 


(^)a^+1>+  Bpl 


+ 


Pp(cos  0)pR1 


(p+2) 


(D. 19c) 


Ap2  '  [T]  Ap3  +  M  Bp2 


(  D .  1 9d  ) 


A 


P1 


=  A 


P2 


B 

p2  1 


(D . 19e) 


REDUCTION  OF  THE  MAGNETIC  POTENTIAL  WHEN  y2  EQUALS  hl  AND  J1(0)  =  0 

The  coefficients  A  3,  Ap2,  Ap3,  Ap4 ,  Bp2,  Bp3>  Bp4,  and  Bp5  for  the  potentials 
are  now  evaluated  for  the  system  consisting  of  a  ferromagnetic  shell  with 
permeability  p2  surrounded  by  an  infinitesimally  thin  current  band  (J2)  in  a  homo¬ 
geneous  medium  with  permeability  in  the  limit  as  p9  =  =  0).  The  varia¬ 

bles  are  defined  in  Figure  9  located  in  the  text.  When  p2  is  set  equal  to  the 
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J 


problem  reduces  to  that  of  finding  the  potentials  in  the  two  regions  of  a  simple 
current  band  (see  Appendix  B  of  Reference  7),  because  the  ferromagnetic  shell  will 
now  have  a  permeability  equal  to  that  of  the  homogeneous  medium  with  permeability 

yl* 

In  this  limit  the  coefficients  should  assume  the  following  form: 


(D.20a) 


3  ,  =  B 
p3  P* 


(D.20b) 


V  (r4"<2p+1)  ) 


(D.20c) 


and  where  A  .  and  B  r  should  reduce  to  coefficients  for  the  potentials  in  the  two 
pi  P5  ? 

regions  for  the  spherical  band  problem.  If  the  coefficients  assume  this  mathemati¬ 
cal  form  it  will  prove  that  the  mathematical  form  of  the  coefficients  for  the  spheri¬ 
cal  shell  surrounded  by  a  thin  current  band  are  mathematically  correct. 

From  Equation  (D.10)  B ^  is 


RP 

R2 

H  M 

rP 

R2 

H  W  [T1 

- 

M  M 

rP  _ 
R2 

[T 

r-(P+D 

R2 

■si 

>]  tTJ 

- 

R2_ 

(p+1) 

rP 

R2 

>] 

'  (D.21) 


Now 


p3 


y2=yl 


Jlimit  B  ,  =  0 

p3 

y2=yi 


(D.22) 
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because 


I^Yj  =  0  for  J 


0 


limit  [T]-  1 


U2=Vil 


The  expression  for  B  .  is  zero  because  J  ,(9)  =  0.  The  expression  for  A 

pZ  P-1-  ] 


(see  Equation  (D.5))  is 


V  ‘  V  M +  Bp2  W 


because 


=  i  and  B 


P2 


m2=m 


1 


0 


C 


The  expression  for  A  .. 

pi 


(see  Equation  (D.18)) 


is 


P1 


=  A 


B  0R 

p2 


-(2p+l) 


and  because  B  _  =  0 


The  expression  for  A 


P3 


(see  Equation  (D.7))  is 


Ap3  -  v  ["]  +  MO] 


(D.27) 


and  because 


=  0 


P2=U 


1 


(D.28) 


Now 


=  A 


=  A  _  =  A 


pi  P2  p3  p4 


^J_£2_(9) 


(2p+l)R^P_1)P1(cos  0) 


(D.29) 


The  expression  for  B 


P4 


(see  Equation  (D.16)) 


is 


3p4 


=  A  V2p+1> 
P3  3 


+  B 


P3  “  p4 


tC2p+l) 

3 


(D.30) 


Because  in  the  limit  = 


and  because  B 


P3 


0  we  have 


(D.31) 
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The  expression  for  B^,.  (see  Equation  (D.17))  is 


B  ,  =  -  A  , 

p5  p4 


(2p+l)  +  B  +  MlJp2(6)  .  .  .... 

4  p4  pR;(p+2)pJ(cos  e) 


(D.32) 


In  the  limit 


B  ,  =  £imit  B  ^ 

p5  p5 


y2~Pl 


VM1 


Vp2(6) 

R“^p+2\2p+l)Pp(cos  6) 


(D.33) 


This  means  that  in  the  five  regions,  the  potential  used  in  Equations  (D.l)  reduce, 
when  =  and  Jpl  =  0,  to  the  form 


W,  II,  III,  IV 


A  rP  ^P^Ccos  0) 

P1  P 


(D. 34a) 


=  I  B 


“TT  \  P1(cos  9) 
-P+1  P 


(D .34b) 


because 


V  -  v  =  V  -  V  =  (2p+1)R(pS)pl(cos  G) 


(D. 35a) 
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m 


p2 


=  B 


p3  p* 


=  0 


(D.35b) 


^3^(0) _ 

P5  R~^P+2)(2p+l)p^  (cos  0) 


(D.35c) 


For  comparison,  the  coefficients  in  Reference  7  are  (primes  are  used  for  distinc¬ 
tion) 


BP2R 


(-2p-l) 

1 


(D.36a) 


B 


I 

p2 


(6) 

_ ±_£ _ 

R,  ^p+^(2p+l)P^(cos  0) 
1  P 


(D.36b) 


When  making  the  comparison  one  lets  R^  =  R^. 

REDUCTION  OF  THE  MAGNETIC  POTENTIAL  WHEN  y2  EQUALS  AND  J2(0)  =  0 

In  a  manner  similar  to  that  of  the  preceeding  section,  the  coefficients  are 
evaluated  for  the  system  consisting  of  an  infinitesimally  thin  current  band[j^(0)] 
surrounded  by  a  ferromagnetic  shell  with  permeability  y2,  in  a  homogeneous  medium 
with  permeability  y^,  in  the  limit  as  y2  =  y^. 

In  this  limit  the  coefficients  should  assume  the  following  form 


=  A  „  =  A 


p2  P3 


p4 


B 

p2  1 


=  B 


=  B 


p2  p3  p4  p5 


(D.37a) 


(D.37b) 


(D.37c) 
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From  Equation  (D.5),  A  is 

P2 


P2 


Vr 


p+1)+(^-vv 


i  + 


,Up+1) 

2 


P2 


P2=|1i 


=  (limit  A  =  A 

p2  p3 

,j2=M1 


The  expression  for  B  „  (see  Equation  (D.6))  is 

P2 


“iV8* 


p2  (2p+l)  (cos  0) 


The  expression  for  A  ^  (see  Equation  (D.4))  is 


M1je2 


(2p+l)R^ 


(  P — 1 ) 


(0) 

P^(cos  0) 
p 


and  because  .1  -  i) 

rT  / 


A 


p4 


0 


Tli 


I’XJU  i 


ion  !  or  B 


Pi 


(see  Equation  (D.10))  reduces  when  p  =  p.  to 


(D. 38) 


(D.39) 


(D .40) 


(D.41) 


(D.42) 
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P3 


(2P+dr: 


Vpl(9) 


P  (cos  0) 
P 


m2=ui 


because,  when  y ^  =  y^and  Jp2  =  0 


=  0 


(D.43) 


(D.44a) 


V_eL(6) 


(2p+l)R~<P+2^P1(cos  0) 


(D ,44b) 


M-1 

(D.44c) 

M  • 0 

(D.44d) 

M  - 0 

(D.44e) 

The  expression  for  A 


p3 


(Equation  (D.7)) 


is 


Ap3  ■  v  [“] + H[>: 


(D.45) 


and  because  f"wl  =  ("xl  =  0  when  =  U1  ,  and  J  „  =  0,  we  have 


Similarly, 


for 


V  '  Ap3  M  +  V  [*] 


(D.47) 


and  because  A  -  =  0  and 

p3 


M 


0  when  ^ 1  *  we  have 


(D.48) 


The  expression  for 


(Equation  (D.18)) 


is 


+  B 


p2 


R- (2p+l) 


(D.49) 


Because  A  =  0,  Equation  (D.49)  yields 
P2 


B  ?R 

p2 


~(2p+l) 

1 


(D.50) 


For  B  ,  , 
p4 


using  Equations  (D.16), 


(D.46) ,  and  (D.42)  we  have 


3  ,  =  B  9 
p4  p3 


(D.51) 


Similarly,  frr  B  ^ ,  using  Equations  (D.17)  and  (D.42),  and  the  fact  that  =  0, 
we  have 
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(D.52) 


B 


P5 


Thus,  we  have  shown  that  when  p. 


yi  and  Jp2(0) 


0, 


A  „  =  A  .  =  A  , 
p2  p3  p4 


0 


(D.53a) 


B  -(2y+1) 

p2  1 


(D.53b) 


B 

P 


2 


B 


P3  '  Bp4  ~  Bp5 


yiJpl<6> 


(2p+l)R~(p+2)P1(cos  0) 


(D ,53c) 


Once  again  this  means  that,  in  the  five  regions,  the  potentials  used  in  Equations 

(D,l)  reduce  when  p,  =  p„  and  J  =  0  to  the  form 
1  l  p2 


v-  £(v  'p)pJ< 


(cos  0) 


(D.54a) 


s 


II, III, IV,  V 


)pp(cos  0) 


y2=Ul 


(D ,54b) 
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For  comparison,  the  coefficients  in  Reference  7  (primes  are  used  for  distinction) 


a',  -  b'  r<-2>- « 

pi  p2  1 


(D.55a) 


MfJ  (0) 


R-(p+2)(2p+i)p1(cos  0) 


(D.55b) 


APPENDIX  E 


FERROMAGNETIC  PROLATE  SPHEROIDAL  BODIES  IN  A  CONSTANT 
EXTERNAL  INDUCING  FIELD 

INTRODUCTION 

8  9 

In  previous  work  Brown  and  Baker  ’  derived  the  closed  form  mathematical 
expressions  for  the  magnetic  flux  density  for  various  configurations  of  a  ferromag¬ 
netic  spheroidal  body  surrounding  and/or  surrounded  by  a  stationary  current  band  of 
azimuthal  symmetry.  The  problem  of  determining  the  magnetic  induction  for  a  prolate 
spheroidal  body  surrounding  and/or  surrounded  by  an  infinitesimally  thin  current 
band  can  be  generalized  to  include  an  external  magnetic  field.  The  superposition 
principle  discussed  in  the  text  of  this  report  can  be  used  in  these  cases  to  include 
a  constant  external  magnetic  field.  The  magnetic  induction  in  each  region  for  a 
three-dimensional  magnetic  spheroidal  shell  in  an  arbitrary  external  magnetic  field 
is  added  to  the  magnetic  induction  for  the  corresponding  region  for  the  spheroidal 
shell  surrounding  and/or  surrounded  by  a  stationary  current  band.  The  problem  of 
deriving  the  magnetic  induction  for  a  current  band  of  finite  width  surrounding  a 
solid  ferromagnetic  spheroid  can  also  be  generalized  to  include  an  external  magnetic 
field  Hq  in  a  similar  manner.  Thus,  the  magnetic  induction  for  a  ferromagnetic 
spheroidal  body  in  an  external  magnetic  field  must  be  determined. 

Both  constant  external  field  problems  were  solved  by  Nixon^  of  the  Center.  The 
closed  form  mathematical  solutions  for  the  magnetic  induction  for  both  constant 
external  field  problems  were  presented  in  Reference  6  in  Cartesian  coordinates.  It 
was  necessary  to  convert  these  mathematical  expressions  to  spheroidal  coordinates 
to  be  compatible  with  this  work. 

SOLID  FERROMAGNETIC  SPHEROID  IN  AN  EXTERNAL  INDUCING  FIELD 

The  solid  ferromagnetic  prolate  spheroid  in  a  constant  external  inducing  field 
is  shown  in  Figur’’  (E.l).  The  permeability  of  the  solid  spheroid  is  and  the 
boundary  of  the  is  determined  by  n  =  =  constant.  The  permeability  p^  of 

vacuum  that  is  external  to  the  spheroid  is  denoted  by  p^.  The  constant  arbitrary 
magnetic  field  is  designated  as  HQ» 

Lt  is  assumed  that  in  the  spheroid  is  constant,  and  that  is  constant  in 

the  region  external  to  the  spheroid.  Because  there  are  no  currents  in  any  regions 
in  the  problem,  the  magnetic  field  H  can  be  expressed  as  the  negative  of  the 
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gradient  of  a  magnetic  scalar  potential  <j>  in  regions  I  and  II,  respectively. 
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(E.ld) 


The  major  step  toward  solving  this  problem  is  to  determine  the  solutions  of  the 
scalar  Laplace's  equation  in  regions  I  and  II  which  satisfy  the  boundary  conditions 
at  n  =  rn  •  In  terms  of  B  and  H,  the  magnetostatic  boundary  conditions  are 


‘hi-V- 


12 


=  0 


at  n  =  n. 


(E .2a) 


ni2 x  (HirV  =  °  at  T\  =  n1 


(E. 2b) 


where  n^0  is  the  unit  vector  normal  to  the  surface  of  the  spheroid,  outward  from 
region  I  to  region  II.  The  general  expression  of  <J>  and  which  satisfy 

Laplace's  equation  in  regions  I  and  II,  are: 
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where  a  is  one-half  of  the  focal  length. 

The  coefficients  are  determined  by  the  magnetostatic  boundary  conditions 
Equations  (E,2a)  and  (E.2b). 
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For  details  of  this  derivation  t  tie  reader  should  consult  Nixon.  The  authors  have 
changed  his  mathematical  expressions  for  the  magnetic  induction  in  regions  I  and  II 
from  Cartesian  to  spheroidal  coordinates.  This  makes  the  expressions  for  the  mag¬ 
netic  induction  compatible  with  the  work  presented  in  the  text  of  this  report. 

These  mathematical  expressions  are,  in  regions  1  and  II, 
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The  problem  of  the  sphe'roidal  shell  is  similar  to  the  problem  of  the  solid 
spheroid.  The  inner  boundary  of  the  prolate  spheroidal  is  and  the  outer  boundary 
is  n2  (see  Figure  E.2)).  The  permeability  of  the  magnetic  material  in  the  shell  is 
U2  and  the  permeability  of  vacuum  that  is  internal  and  external  to  the  shell  is 
denoted  by  p. .  The  constant  external  magnetic  field  is  designated  by  Hq. 

The  problem  of  deriving  the  closed  form  mathematical  expressions  for  the  mag¬ 
netic  flux  density  in  each  of  the  three  regions  (I  through  III)  was  solved  in  detail 
by  Nixon. ^  The  problem  was  solved  in  a  method  exactly  analogous  to  the  method  used 
to  solve  the  problem  of  a  solid  spheroid  in  a  constant  external  magnetic  field. 

For  details  of  the  derivation,  Reference  6  should  be  consulted. 

The  general  expressions  for  the  magnetic  scalar  potential  <p  in  regions  I 

m 

through  III  are 
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The  coefficients  are  determined  by  the  usual  magnetostatic  boundary  conditions 
on  the  spheroidal  surfaces  at  r]  =  n1  and  n  =  n2  •  The  coefficients  determined  by 
this  method  are 


(E.7) 
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(Note:  Above  equation  continued  on  next  pa<  .) . 
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(Note:  Above  equation  continued  on  next  page). 
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Using  the  symmetry  conditions  that  exist  in  this  problem,  the  last  four 
stants  (M  ,  F  ,  K  ,  and  L  )  are  obtained  from  the  previous  equations  by  simple 
substitution.  Therefore,  M  ,  F  ,  K  ,  and  \!  are  determined  by  substituting  H 
Hox  in  Equations  (E.12)  through  (E.16). 


The  authors  have  changed  Nixon’s  mathematical  expressions  for  the  magnetic 
induction  in  regions  I  through  III  from  Cartesian  to  spheroidal  coordinates.  This 
makes  the  expressions  for  the  magnetic  induction  compatible  with  the  work  presented 
in  the  text  of  this  report.  The  mathematical  expressions  in  regions  I  through  III 
are 
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(Note:  Above  equation  is  continued  on  next  page.) 
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